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Pitch  angle;  cot  0  =  2Tra/p 


l  ,n 

space 

space 


xi  i  i 


Wavelength 


\ 


S 

A 

♦ 

♦ 

a) 


A  series  summation  defined  in  Eq.  E.3 
Permeability;  integer  index 
Permeability  in  Region  i 
Permeability  of  free  space 
Integer  index 
Integer  index 

Angular  tape  width,  £  -  2ir5/p 
Unit  vector  in  the  $  direction 
Dependent  variable 
Angular  frequency 


*  Cross  product  of  two  vectors 

*  Dot  product  of  two  vectors 

*  Sign  denoting  multiplication;  superscript  denoting  com¬ 
plex  conjugate 


I.  INTRODUCTION 


1 .1  History 

With  the  invention  of  the  traveling-wave  tube  (TWT)  by  Rudolf 
Korapfner  [1]  in  19^2,  the  helix  structure  was  established  as  the  delay 
mechanism  used  to  slow  down  the  forward  propagation  of  the  electromag¬ 
netic  waves.  It  was  this  same  helical  geometry  with  its  diner  jr.tial 
screw  symmetry  which  ultimately  proved  to  yield  the  least  dispersion  and 
widest  bandwidth.  However,  a  problem  was  encountered  for  high  voltage 
operation:  backward-wave  oscillation  (BWO).  At  high  voltages,  the 
axial  focusing  of  the  electron  beam  became  difficult,  which  resulted  in 
the  beam  interacting  with  field  components  other  than  the  fundamental. 
At  a  certain  voltage,  interaction  with  the  space  harmonics  became  strong 
enough  to  induce  backward-wave  oscillation.  To  compound  the  problem, 
the  impedance  for  the  electron  interaction  with  the  fundamental  com¬ 
ponent  of  the  fields  was  reduced  because  of  the  increased  energy  content 
of  the  noninteracting  space  harmonics.  Though  advanced  focusing  methods 
helped  to  resolve  these  problems,  a  better  solution  was  a  device  which 
had  a  larger  interaction  impedance  of  the  fundamental  component  relative 
to  the  space  harmonics,  while  maintaining  the  wide  band  characteristics 
of  the  helix.  The  contrawound  helix  proved  to  be  such  a  device. 

The  contrawound  helix  shown  in  Fig.  1.1a  was  first  investigated  by 
Chodorow  and  Chu  [2]  in  195^.  They  observed  that  such  a  structure, 
consisting  of  two  tape  helices  wound  in  opposite  directions,  could  be 
qualitatively  analyzed  by  considering  the  simple  superposition  of  the 
two  single-helix  fields.  In  one  situation,  the  fields  were  thought  of 
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Fig.  1.1.  The  two  tape  contrawound  helix  (a)  and  the 

ring-bar  version  of  the  oontrawound  helix  (b). 


as  being  superimposed  180°  out  of  phase.  This  yielded  a  field  configu¬ 
ration,  labeled  the  antisymmetric  mode,  with  a  decreased  axial  field  and 
an  increased  radial  field  relative  to  the  usual  single  helix  fields.  A 
second  mode,  the  symmetric  mode,  resulted  when  the  two  single-helix 
modes  were  considered  superimposed  in  phase.  For  such  a  mode,  the  axial 
electric  fields  of  the  fundamental  component  added  rather  than  sub¬ 
tracted,  giving  a  stronger  axial  field,  while  at  the  same  time  reducing 
the  radial  field.  Furthermore,  the  stored  energy  associated  with  the 
fundamental  component  of  the  magnetic  field  —  energy  which  is  useless 
for  electron  beam  interaction  —  was  found  to  approach  rsro.  The  impli¬ 
cation  was  that  the  fundamental  component  of  the  symmetric  mode  primarly 
carried  electric  energy  and  that  the  space  harmonics  carried  principally 
magnetic  energy.  This  was  in  contrast  to  the  single  helix  in  which  the 
electric  and  magnetic  energy  were  roughly  equal  in  the  fundamental 
component  of  the  operating  mode.  However,  such  qualitative  analysis 
ignored  completely  the  interaction  between  the  two  helices. 

The  quantitative  analysis  of  the  contrawound  helix  performed  by 
Chodorow  and  Chu  was  simply  that  of  solving  Maxwell’s  equations  under 
the  appropriate  boundary  conditions,  i.e.,  a  boundary  value  problem. 

i 

The  difficulty  arose  in  defining  adequately  these  boundary  conditions 
,  and  describing  correctly  the  field  configuration  for  the  particular  mode 

under  investigation.  To  accomplish  this,  the  electric  and  magnetic 
|  fields  were  expressed  in  terms  of  the  surface  current  density  on  the 

i 

helices.  Together  with  a  technique  involving  variational  calculus, 
these  expressions  were  used  to  obtain  a  determinantal  equation  in  which 
the  propagation  constant  in  free  space  was  written  as  a  function  of  the 


-  3  - 


axial  propagation  constant  of  the  fundamental  space  harmonic.  The 
variational  method  was  chosen  because  it  had  the  great  advantage  that  it 
permitted  the  use  of  successive  approximations,  which  in  the  final  form 
converge  to  the  exact  solution.  The  result  was  a  determlnantal  equation 
which  was  both  compact  and  numerically  economical. 

Building  on  the  work  of  Chodorow  and  Chu,  Ayers  and  Kirstein  C 3 ] 
examined  the  ring-bar  circuit,  an  easy-to-make  version  of  the  contra- 
wound  helix.  The  ring-bar  structure  illustrated  in  Fig.  1.1b  consisted 
of  a  series  of  rings  connected  one  to  another  by  bars  at  alternate  ends 
of  a  diameter.  To  manufacture,  it  was  simply  a  matter  of  making  a 
number  of  saw  cuts  in  a  tube.  Numerical  analysis  was  done  based  on  an 
unpublished  determinantal  equation  for  the  ring-bar  structure  derived  by 
Chodorow  and  Chu.  This  determinantal  equation  was  completely  unlike 
that  for  the  contrawound  helix.  Rather  than  describing  the  fields  in 
terms  of  the  surface  currents,  assumptions  were  made  about  the  form  of 
the  fields  themselves.  These  were  then  manipulated  by  a  variational 
technique  into  a  determinantal  equation.  Furthermore,  Floquet's  theorem 
was  applied  to  the  step  screw  symmetry  of  the  problem  so  that  mathemati¬ 
cal  orthogonality  was  defined  over  only  half  of  the  normal  period, 
rather  than  the  usual  full  period."^  Useful  results  were  obtained  using 
this  model,  which  matched  experimental  data  quite  well  over  the  first 
portion  of  the  dispersion  diagram. 


t 


Refer  to  Appendix  B  for  more  details. 


Ayers  and  Kirstein’s  work  was  not  strictly  numerical,  however. 
Mostly,  they  dealt  with  the  experimental  aspect  of  the  problem  as  did 


others  at  that  time. 

Concurrently,  similar  experimental  work  was  being  performed  on  the 
contrawound  helix  as  well  as  its  related  circuits.  Birdsall  and 
Everhard  [4]  analyzed  various  forms  of  these  circuits  and  how  they  were 
affected  by  such  things  as  dielectric  loading,  helix-to-waveguide  tran¬ 
sitions,  and  periodic  support  stubs.  Nevins  [5]  considered  the  effects 
of  altering  various  geometric  parameters  along  with  examining  the  elec¬ 
tron  beam  interaction.  However,  after  the  initial  flurry  of  work  in 
this  area,  interest  in  contrawound  helix  structures  declined. 
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1  .2  Purpose  of  the  Report 


In  spite  of  their  many  excellent  properties,  contrawound  helix 
type  circuits  have  been  neglected  in  the  years  since  their  initial 
development.  This  is  mainly  because  they  have  been  difficult  to  manu¬ 
facture  to  the  high  tolerances  necessary,  but  also  because  single-helix 
technology  is  well  established.  However,  with  today's  manufacturing 
capability,  renewed  investigation  into  these  structures  is  warranted. 

It  is  the  purpose  of  this  report  to  reanalyze  the  work  first 
carried  out  by  Chodorow  and  Chu  [2]  on  an  unloaded  contrawound  helix. 
This  mathematical  model  is  then  extended  to  include  both  the  effects  of 
a  surrounding  conducting  sheath  as  well  as  dielectric  loading.  The 
results  obtained  are  then  compared  to  previously  published  experimental 
results . 
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1 . 3  Organization 


This  report  deals  primarily  with  Chodorow  and  Chu's  [2]  contra- 
wound  helix  circuit  in  free  space  and  the  more  general  problem  in  which 
shielding  and  dielectric  loading  are  considered. 

Chapter  Two  investigates  the  free  space  problem.  It  begins  with  a 
discussion  of  the  mathematical  formulation  used  to  derive  the  dispersion 
equation.  It  then  proceeds  to  comment  on  the  numerical  results  and  how 
they  compare  to  experimental  results  for  similar  cases. 

Chapter  Three  is  concerned  with  the  general  problem  of  a  contra- 
wound  helix  surrounded  by  a  conducting  sheath  and  dielectrically 
loaded.  It  is  shown  how  these  new  boundary  conditions  effectively  alter 
the  dispersion  equation  and  how  this  alteration  affects  the  numerical 
solution.  Furthermore,  these  results  are  compared  with  the  appropriate 
experimental  results  so  as  to  determine  the  effectiveness  of  the  changes 
to  the  determinantal  equation. 

Finally,  several  topics  related  to  the  solution  of  the  contrawound 
helix  boundary  value  problem  are  covered  in  the  appendices:  the  general 
form  of  Floquet's  theorem  for  step-screw  periodicity,  the  general 
Fourier  expansion  of  the  electromagnetic  field  functions  in  circular 
cylindrical  coordinates,  and  a  discussion  of  the  associated  Lagrangian 
which  Is  used  in  conjunction  with  a  variational  method  to  obtain  a  solu¬ 
tion. 
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1.4  Mathematical  Preliminaries 


The  mathematical  problem  is  simply  to  solve  Maxwell's  equations 
under  the  appropriate  boundary  conditions,  i.e.,  a  boundary  value  prob¬ 
lem.  The  problem  is  constrained  to  be  source  free,  so  that  for  the 
electric  field  (E) , 

V  •  E  -  0  (1.1) 

everywhere  in  space.  Further  conditions  are  that  the  structure  is 
lossless  and  that  the  solutions  are  restricted  to  the  time  harmcr.i; 
form  e  These  restrictions  allow  Maxwell’s  equations  to  be  manipu¬ 

lated  into  wave  equations  for  the  electric  field, 

V2E  *  k2E  -  0  (1.2) 

and  for  the  magnetic  field  (H), 

V2H  +  k2H  =  0  (1.3) 

where 

k  -  uj/yT  ( "  .  - ) 

in  which  y  and  e  are  the  permeability  and  permittivity,  respectively,  of 
the  medium  under  consideration,  and  u>  is  the  radian  frequency.  Though 
Eqs.  1.2  and  1.3  are  vector  equations,  they  reduce  to  the  scalar 
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Helmholtz  equation  as  one  considers  separately  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  polarizations  of  the  fields.  However, 
unlike  the  traditional  waveguide  in  which  the  boundary  conditions  can  be 
fulfilled  by  either  a  transverse  electric  or  a  transverse  magnetic 
field,  the  geometry  of  any  slow-wave  structure  is  such  that  a  superposi¬ 
tion  of  these  two  fields  is  necessary  to  satisfy  the  boundary  condi¬ 
tions  . 

1.4.1  Symmetry  and  Periodicity 

The  contrawound  helix  has  several  symmetry  characteristics 

uniquely  associated  with  this  class  of  structures.  There  are  two  planes 

of  reflective  symmetry,  the  (r,  8}  plane  and  the  (r,  z)  plane,  each 

intersecting  at  the  crossover  point  of  the  two  helices.  Considered 

together  with  Maxwell's  equations,  these  reflection  symmetries  require 

that  any  solution  must  be  either  even  or  odd  in  z  and  <j> ,  excluding 

degeneracy.  Such  solutions  are  standing  waves  which  can  be  combined  to 

give  running  waves.  Each  plane  of  reflective  symmetry  yields  two  types 

of  solutions,  depending  on  whether  a  conducting  or  magnetic  wall  is 

considered.  The  result  is  a  total  of  f our  types  of  field  configurations 

(Table  l.i).  Of  interest  is  the  scena-io  in  which  the  (r,  z)  plane  of 

reflective  symmetry  is  replaced  by  a  magnetic  wall.  The  result  is  that 

the  vector  component  of  the  electric  field  in  the  z  direction,  E  ,  is 

-z 

even  in  $,  while  that  for  the  magnetic  field,  H  ,  is  odd  in  <j>.  Because 
Hz  is  odd  in  $ ,  it  must  vanish  on  the  axis.  In  particular,  the  funda¬ 
mental  Fourier  component  of  is  identically  zero  and  the  energy  asso¬ 
ciated  with  the  fundamental  space  harmonic  for  the  system  is  stored 


Table  1.*.  Symmetry  of  the  field  components  for  the 
the  symmetric  and  antisymmetric  modes. 


principally  In  the  electric  field.  The  combination  of  solution  types  1 
and  2  produces  a  wave  traveling  in  the  z  direction  with  these  desirable 
properties  and  is  called  the  symmetric  mode. 

The  antisymmetric  mode  is  a  combination  of  the  basic  field  types  3 
and  This  yields  another  traveling  wave,  but  one  in  which  the  elec¬ 
tric  field  is  shorted  out  along  the  z  axis.  Because  the  (r,  z)  plane  of 
symmetry  is  now  a  conducting  wall,  the  energy  associated  with  the  fun¬ 
damental  component  is  stored  primarily  in  the  magnetic  field,  a  char¬ 
acteristic  detrimental  for  TWT  operation. 

Like  all  periodic  structures,  the  contrawound  helix  is  invariant 
under  the  transformation , 

(r,  4i,  z)  (r ,  <|>,  z  ±  p)  (1.5) 

where  p  is  the  period.  This  relationship  requires  a  certain  functional 
dependence  in  the  description  of  the  fields  for  the  structure.  Known  as 
Floquet’s  theorem,  this  requirement  states  that  under  a  translation  of 
an  integral  number  of  periods,  the  fields  can  differ  at  most  by  a  con¬ 
stant.  The  result  is  that  the  axial  propagation  characteristics  are 
limited  to  a  particular  form. 

Floquet's  theorem  is  also  applicable  to  step-turn  periodicity, 
also  known  as  screw  symmetry,1’  For  the  case  of  the  single  helix  (Fig. 
1.2),  the  step-turn  periodicity  is  described  by  the  differential  screw 

*  See  Appendix  B. 
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transformation, 


(r,  *,  z) 


z  ♦  6z 


) 


(1.6) 


where  the  choice  of  sign  depends  on  whether  the  helix  is  right-  or  left- 
handed.  The  form  of  the  fields  remains  invariant,  again  to  within  a 
constant,  for  a  differential  step  in  the  z  direction  coupled  with  the 
appropriate  amount  of  differential  twist.  The  situation  with  the  con- 
trawound  helix  is  that  this  "differential"  step  is  uniquely  half  of  the 
axial  period,  6z  -  p/2,  such  that  transformation  of  Eq.  1.6  becomes 

r,(j>±Tt,z+|'j  (1.7) 


(r,  $,  z) 


1 

! 

I 

I 


As  with  the  traditional  application  of  Floquet's  theorem,  the  step-turn 
symmetry  imposes  further  restrictions  on  the  propagation  characteris¬ 
tics.  As  one  might  expect,  these  particular  restrictions  are  an  aspect 
of  coupling  between  the  z  and  $  coordinates.^ 


1.^.2  Boundary  Conditions  and  Space  Harmonics 

It  is  no  simple  matter  to  satisfy  the  boundary  conditions  for  the 
contrawound  helix.  A  comparison  with  the  circular  cylindrical  waveguide 
shows  the  inherent  difficulties  of  matching  the  boundary  conditions  for 
a  slow  wave  structure.  The  functional  form  of  the  field  intensities  for 


*  See  Appendix  B. 


a  perfectly  conducting  waveguide  consists  of  regular  cylindrical 

I 

I 

Bessel's  functions,  each  having  an  infinite  number  of  zeros.  This 
property  enables  each  mode  to  individually  satisfy  the  boundary  condi¬ 
tions  and  to  therefore  exist  independently  of  all  the  other  modes.  It 
is  the  closed  nature  of  the  waveguide  boundary  which  allows  this  type  of 

i 

separation  among  the  modes.  1 

In  the  mathematical  description  of  the  field  configurations  for 
the  slow  wave  structure,  the  regular  Bessel's  functions  are  replaced  by  1 

modified  Bessel's  functions  which  have  no  zeros.  Thus,  it  is  not 
possible  to  satisfy  the  boundary  conditions  uniquely  for  each  mode,  but 
rather  the  solution  is  found  in  an  aggregate  of  these  modes.  The  dis¬ 
tinction  is  then  made  that  these  "modes"  are  not  really  modes  in  the  i 

sense  that  they  can  exist  independently  of  each  other,  but  are  instead 

waves,  termed  space  harmonics,  which  must  exist  in  unison  to  satisfy  the  1 

I 

boundary  conditions.  These  space  harmonics  are  related  by  the  period- 

< 

icity  of  the  structure,  and  each  is  orthogonal  to  the  rest,  in  r,  <p ,  and 
z.  For  the  contrawound  helix,  the  $  dependence  of  each  harmonic  either 
has  the  form  cos  (n<t> )  or  sin  (n$).  The  component  with  n  *  0  is  labeled 
the  fundamental  space  harmonic. 

1.4.3  Dispersion  Equations  and  Solution 
by  Variational  Calculus 

With  Maxwell’s  equations  satisfied  and  the  boundary  conditions 
correctly  accounted  for,  a  determinants!  equation  can  be  found  for  the  I 


slow  wave  structure.  The  dispersive  characteristi cs  —  how  frequency 
varies  as  a  function  of  the  phase  constant  (phase  velocity  over 


frequency)  —  are  described  by  this  equation.  It  is  obtained  formally 
by  the  usual  technique  of  analysis  in  terms  of  orthogonal  functions 
which  yield  four  infinite  sets  of  homogeneous  simultaneous  equations. 
The  overall  system  is  then  solved  by  well  known  matrix  methods. 

This  formal  method  of  solution,  however,  is  especially  inconveni¬ 
ent  in  the  case  of  the  contrawound  helix.  For  unlike  the  helix,  the 
simultaneous  equations  describing  the  contrawound  helix  are  doubly 
infinite  over  two  indices.  Thus,  a  variational  technique  is  used  for 
deriving  approximate  solutions  to  yield  numerical  results.  Though  this 
technique  also  leads  to  the  same  infinite  set  of  equations  just 
described,  its  advantage  is  that  it  allows  one  to  systematically  approx¬ 
imate  the  solution  to  the  eigenvalue  problem. 
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II.  THE  CONTRAWOUND  HELIX  IN  FREE  SPACE 


2.1  The  Boundary  Value  Problem 

The  analysis  of  the  contrawound  helix  in  free  space  is  initiated 
by  separating  the  problem  into  two  regions,  one  inside  and  one  outside 
the  cylindrical  surface  r  -  a,  Region  1  and  Region  2,  respectively,  as 
shown  in  Fig.  2.1.  The  helices  are  assumed  to  be  infinitely  long,  of 
equal  radii  (r  -  a),  and  wound  with  an  infinitely  thin  perfectly  con¬ 
ducting  tape.  The  dimensional  quantities  that  describe  the  structure 
are  the  pitch,  the  tape  width,  and  the  radius.  These  quantities  are 
related  as  follows: 


2ira 

P 


cot  e 


2  tt  6 
P 


5 


where 

a  «  helix  radius 
p  -  period  of  helix 
6  =  tape  width 
9  »  pitch  angle 


(2.1) 


(2.2) 
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2.1.1  The  Field  Functions 


To  satisfy  Maxwell's  equations  in  the  form  of  the  Helmholtz  equa¬ 
tion  as  well  as  the  complicated  boundary  conditions,  it  is  necessary  to 
construct  a  solution  using  the  usual  technique  of  Fourier  decomposi¬ 
tion.  Thus,  each  field  component  is  written  as  an  infinite  sum  of 
elementary  waves  in  which  each  wave  satisfies  the  differential  equation 
and  some  of  the  symmetry  conditions.  By  superimposing  the  total  set  of 
waves,  all  other  restricting  conditions  can  be  satisfied. 

Any  arbitarary  field  in  a  homogeneous  source-free  region  can  be 
expressed  as  the  sum  of  a  TM  field  and  a  TE  field.  The  skew  boundary 
conditions  of  the  contrawound  helix  make  it  necessary  to  have  both  TE 
and  TM  fields  present  in  any  given  mode.  Thus,  the  Fourier  decomposi¬ 
tions  describing  Ez  and  Hz  in  the  two  regions  are 
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From  Maxwell’s  equations,  the  other  field  components  are  written  in  each 
region  in  terms  of  E  and  H  These  are  in  Region  1, 


(0  <  r  <  a) 
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and  in  Region  2 
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In  Eqs .  2.3  through  2.10,  I^^r)  and  K£(Ynr)  are  modified  cylindrical 
Bessel  functions.  These  are  necessary  for  slow  waves  (phase  velocity 


Refer  to  Appendix  A. 
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less  than  c,  the  speed  of  light)  and  are  chosen  to  give  nonradiating 
solutions.  Any  derivative  of  these  functions  is  with  respect  to  r.  The 
relationship  between  Bn  and  is  provided  by  Maxwell's  equations  (again 
for  slow  waves)  and  is 


y 


2 

n 


(2.11  ) 


where  k  Is  the  wave  number  as  defined  in  Eq.  1.4.  From  periodicity  and 
Floquet's  theorem,  Bn  has  the  property, 


B 


n 


2irn 


(2.12) 


such  that  8Qp  is  the  phase  shift  per  period  along  the  structured 
However,  the  screw-symmetry  or  step-turn  periodicity  provides  an  addi¬ 
tional  relationship  between  the  two  summation  indices  l  and  n,  whereby 
the  axial  propagation  constant  can  be  redefined  to  be* 


a.n 


B  ♦ 

0,0 


U 


*")f 


(2.12’  ) 


B 


0,0 


z 


B 


o 


In  light  of  Eq.  2.12,  Eq.  2.11  is  then  reformed  as 


*  Refer  to  Appendix  B  for  details. 
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(2.13) 
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2 

i,n 


k 
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2.1.2  Boundary  Conditions 

The  boundary  conditions  to  be  satisfied  are  that 

a 

n  *  E  is  continuous  everywhere  on  the 

the  cylindrical  surface  r  «  a  (2.i^a) 


except  on  the  helices  where  it  is  zero. 


n 

X  E  - 

0  on  helices 

(2.Ub) 

and  that 

A 

n 

X  H 

is  continuous 
except  on  the 

on  r  »  a, 
helices 

(2.15) 

where  it  is  proportional  to  the  surface  current  density.  In  the  above 

a 

conditions,  n  is  the  unit  vector  normal  to  the  surface  r  =  a. 


2.1.3  The  Field  Components  in  Terms  of 
the  Surface  Current  Densities 


Separating  condition  2.1^a  into  its  orthogonal  components  yields 

the  equations: 


and 


lEz 


r 

2~z 


r  =  a 


(2.16) 


r  -  a 


(2.17) 
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From  Eq.  2.16,  the  Eqa.  2.3  and  2.4  are  equated.  After  applying  the 
principle  of  orthogonality  over  one  period  on  the  cylindrical  surface  r 
-  a,  the  relationship  between  the  two  sets  of  Fourier  coefficients  is 


obtained: 


l,n 


B 

*.n  1^ 


(2. 18) 


where  I  =  I  (y  a),  etc.  Proceeding  similarly  from  Eq.  2.17,  Eqs.  2.7 
it  t  n 

and  2.9  are  employed  to  yield  the  relationship 


Ci,n  *  Dt,n 


(  c. .  r  , 


Examined  next  is  the  boundary  condition  for  the  H  fields,  Eq. 
2.15.  Like  the  fields,  the  surface  current  density  on  the  helices  can 
be  expanded  in  a  convenient  form  for  algebraic  manipulation.  Its  <j>  and 
z  directed  components  are  decomposed  as 


-  I  e  JB*’nZ 

i ,n  l,n 


(2.20) 


•  I  *JI.n  * 


-JB.  Z 


£,n  ^  J 


l  ,n 


(2.2' ) 


V  -  I  ZJ'  e'JSt’"Z  eJt‘ 

l)n  *•" 


(2.22) 


and 


j 


< 


-  22  - 


e 


(2.23) 


V  -  I 


£  ,n 


Zj» 

t,n 


‘J6l,nZ  gJ£4> 


in  which  the  "+"  and  superscripts  designate  the  left-handed  and 
right-handed  helices,  respectively,  and  in  which  £  and  n  are  allowed  to 
range  from  minus  infinity  to  plus  infinity.  In  terms  of  the  H  field 
components,  the  condition  of  Eq.  2.15  yields  the  boundary  equations: 


H  -  „H  -  *J~  ♦  *J+  ,  r  «  a  (2.24) 

1  z  2  z 

and 

„H  -  H  -  ZJ~  ♦  2J+  ,  r  -  a  (2.25) 

2  $  1  <J> 


by  writing  Eq.  2.24  in  terms  of  the  appropriate  Fourier  expansions  (Eqs. 

2.5,  2.6,  2.2C,  and  2.21),  orthogonality  is  used  to  express  the  two  sets 

of  field  expansion  coefficients  (c„  and  D„  )  in  terms  of  the  current 

£ ,  n  i ,  n 

density  coefficients.  This  gives,  after  making  use  of  the  relationship 
in  Eq.  2.19, 


C 


£  ,n 


£,n) 


Vi 


(2.26) 


and 

D  .  ^  [zjLs  : 

Ul,n  K’I4  -  K£I’ 

Finally,  considerations  of  Eq.  2.25  in  which  the  appropriate  expansions 
(Eqs.  2.8,  2.10,  2.22,  and  2.23)  are  substituted  yields,  after  some 


algebra,  the  result  that 


(2.28) 


and 


l 


where 


l 


(2.29) 


v,  _ 

1 ,n  * 

r,  \ 

i  zj"  ♦  v  i 

£3„ 

i,  n 

(V  +  V 

I  i.n  l,nj 

/ 

r  1 

\h*i  " 

ju)E 

\  Ji,n  Jltn] 

Y2 

i,n 

V  J 

L  J 

(2.30) 


Again  note  that  the  modified  Bessel  functions  have  the  argument  Y  ^a. 

X.  f . « 

The  denominator  in  Eq.  2.30  is  the  wronskian  of  (i  ,  K  )  and  can  be 
replaced  by  -“[/(y^  na ) : 


’  riKi 


y.  a 
i  ,n 


(2.3i) 


Having  written  the  Fourier  coefficients  for  the  fields  in  terms  of  the 
current  density  coefficients,  Eqs .  2.26  through  2.29,  the  field  intensi¬ 
ties  are  at  last  expressed  in  terms  of  the  surface  current  density.  And 
since  all  field  quantities  are  defined  to  within  a  constant,  the  summa¬ 
tions  can  be  multiplied  by  the  quantity  -/e//y  so  as  to  express  the 
field  components  in  the  following  form: 


V  •  *■  z)  ■  jk  J  (Yna)2 

Z  ,n 


z  -  z* 
J«  _  +  J 


18.  a 
l  ,n 


*'n  l’nJ  (y,  a)2 

v  l  ,n  ' 


l  l  ,n  ^2,  ,n 


ItKt 


-  2U  - 


(2.32) 
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Y„  a 
2,  ,n 


*  K£  ^a,nr) 


e'JBZ,nZ  eJl* 


(2.35b) 


J 

In  Eqs.  2.32  and  2.33,  the  modified  Bessel  functions  have  for  their 

arguments  Y  a,  and  throughout  Eqs.  2.32  through  2.35,  the  indices  l 

*> ,  n 

and  n  range  from  minus  infinity  to  infinity. 


2.1.4  The  Surface  Current  Densities  on  the  Helices 

With  the  field  components  written  in  terms  of  and  ZJ^  ,  it 

x  i  n  i  n 

is  important  to  analyze  in  some  detail  the  characteristics  of  the  sur¬ 
face  current  densities. 

From  the  boundary  condition  of  Eq.  2.15,  H  and  H  are  continuous 

0  z 

at  r  *  a,  except  on  the  helices.  To  meet  this  condition,  the  surface 
current  densities  must  be  constrained  as 


J  *  0  off  helix  (-) 


(2.36a) 


J  ■  0  off  helix  (♦)  (2.36b' 


26  - 


In  other  words,  there  can  be  no  current  off  the  helix  tapes.  A  con¬ 
straint  such  as  this  is  somewhat  of  a  novelty  in  boundary  value  prob¬ 
lems.  In  this  case  (as  well  as  other  open  structures),  the  boundary  at 
r  -  a  must  be  explicitly  defined  to  satisfy  the  geometrical  considera¬ 
tions.  Similarly,  the  current  amplitudes  are  chosen  to  satisfy  the 
symmetry  properties  and  the  particular  mode  of  operation. 

To  obtain  the  symmetric  mode  as  defined  in  Section  1.U.1,  it  is 

necessary  for  E  ,  E  ,  and  H  to  be  even  in  *,  while  H  ,  H  ,  and  E  are 

-r  -z  -$  -r  -z 

odd  in  $.  From  Eqs.  2.24  and  2.25,  it  follows  that  the  current  densi¬ 
ties  must  satisfy  the  same  symmetry  conditions  as  the  fields.  Thus, 

J  (a,  <(>  z)  ■  J  (a,  -$,  z)  (2.37a) 

and 

^J+(a,  $,  z)  «  -  (a,  ~<j>,  z)  (2.37b) 

To  facilitate  the  narrow  tape  approximation  in  which  it  is  assumed 
current  flow  is  primarily  in  a  direction  parallel  to  the  helices,  it  is 
advantageous  to  express  ZJ±  and  in  terms  of  components  parallel  and 

perpendicular  to  the  appropriate  helix  tape.  Referring  to  Fig.  2.2  and 

noting  that  each  pair  of  currents  ( ,  lj  +  )  and  (^J  ,  )  is  oriented 

symmetrically  with  respect  to  the  z  axis,  the  following  expressions  may 


be  formed: 


Fig. 


sin  e 


cos  e 


(2.38a) 


ZJ  (a,  z) 

-  l«J~(a,  <*>,  z)  ±  XJ  (a,  *,  z) 

*j~(a,  $,  z)  cos  6  sin  6  (2.38b) 

zJ+(a,  <f>,  z)  sin  0  cos  0  (2.39a) 

-  ±  lj*(a,  <p,  z)  ♦  XJ  (a,  4> >  z) 

*J+(a,  if,  z)  cos  6  sin  8  (2.39b) 

Condition  2.37  then  becomes 

(a,  $,  1 )  m  (a,  ~4>,  z)  (2.410a) 

V(a,  z)  -  V(a,  -4>,  z)  (2.H0b) 

With  the  Fourier  coefficients  for  and  J_J±  defined  in  the  usual 

way, 


'Aa,  4»i  z) 
“^(a,  4>,  z) 


I 

i,n 


I 

Z,n 


eJl* 


(2. a ) 


( 2 .  in  b ) 


the  field  expansions  can  be  written  in  terms  of  these  coefficients. 
However,  by  making  the  assumption  of  narrow  helix  tapes,  the  resulting 
expressions  for  the  fields  are  greatly  simplified.  Mathematically,  this 
assumption  translates  to 
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"Sj  cot  8  <<  and  "*\J  tan  6  << 


from  which  the  set  of  Fourier  coefficients 


i.n 


is  seen  to  dominate.  Thus,  Eqs.  2.38  and  2.39  can  be  reformed  as 


z  — 

J  (a,  4>,  z)  sin  6 

»  (a,4>,z) 

(a,  <fr,  z)  cos  0 


(2.42a) 


(2.42b) 


and 


ZJ+(a,  <f>  ,z) 


*  + 

J  (a,  4>,  z) 


sin  0 


±  ^J+(a,$,z) 


cos  0 


(2.42c) 


( 2 . 4  2d ) 


Finally,  since 


6  *  ft 

,n 


6_,  -  6.  , 

x  f  n  x  f  n  x. 


(2.43a) 


(2.43b) 
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I 


the  symmetrical  relationships  between  the  Fourier  coefficients  for  *j+ 
and  can  be  established.*  Making  use  of  Eqs.  2.M0a  and  2. Mia,  t^ese 
become 


£,n 


lr 


-£,n+£ 


(2. MMa) 
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J-l,n+£ 


(2. MMb) 


2.1.5  The  Determinantal  Equation  from  the  Exact 
Solution  to  the  Boundary  Value  Problem 

As  mentioned  in  the  previous  section,  because  of  the  nonhomogene- 
ous  character  of  the  boundary  at  r  *  a,  the  geometry  of  the  structure 
must  be  reflected  in  the  nature  of  the  surface  current  densities. 
Specifically,  each  surface  current  density,  J+  and  J  ,  cannot  exist  off 
its  respective  helix  (Eq.  2.36).  And  from  the  relation  between  H  and  J, 
it  follows  that  this  is  also  a  restriction  on  the  H  fields  at  r  -  a. 
There  is  a  similar  condition  for  the  E  fields;  namely,  that 


E  (r  -  a)  •  0 
z 

on  helices 

(2.M5a) 

E  (r  -  a)  -  0 

on  helices 

(2.M5b) 

which  is,  in  fact,  the  boundary  condition  2.1Mb. 

To  satisfy  Eq.  2.36,  the  constraint  on  the  surface  current  densi¬ 
ties,  the  usual  technique  of  analysis  in  terms  of  orthonormal  functions 


*  See  Appendix  D  for  details. 
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is  used.*  Needed  are  two  complete  sets  of  orthonormal  functions 

(labeled  g+  and  g  1,  one  for  each  helix  tape,  defined  on  the  tapes 
y,v  °y,v 

Z  +  i  i 

only.  Making  use  of  orthogonality  allows  J"  ana  J.  to  be  written 

*>  f  n  1 1  n 

in  terms  of  these  functions  so  that  each  of  these  current  coefficients 

will  be  compatible  with  the  constraint  of  Eq.  2.36. 

The  result  is  that  the  E_  and  E  expressions  are  now  summed  over 

z  4> 

four  indices  rather  than  just  two,  and  Eqs.  2.^5a  and  2.li5b  take  the 
form 


E  (r  -  a)  - 
z 


i 


4,n;y,v 


^  Jy ,  v  Gt,n;u,v  "  ^Bl,na^ 

y  9  V 


I  V  g7  ♦  l  V  g! 

y,v  t,n;y,v  L  y,v  l,n;y,' 
y  ,  v  ’  y  ,  v 


J  J 


»  P‘J6i.nz  Jl<p 


e Y  -  0  on  helices 


(2.U6a) 
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l 
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Refer  to  Reference  2  for  detailed  discussion,  pp.  38-^6. 
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*  0<a)8  KJIJ 
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-  I  V  ot 

lifv  £ ,  n ;  u  ,  v 

y .  v 


e’JBl.n?  eJU  .  0 


on  helices 


(2.46b) 


In  Eqs.  2.46a  and  2.46b,  the  functions  G*  and  G„  are  the 

£,n;y,v  l ,n; y , v 

result  of  applying  orthogonality  to  the  two  sets  of  orthonormal  func¬ 
tions  g+  and  g~  . 

y,V  UtV 

Finally,  one  can  operate  E_  and  E  with  the  orthonorraal  functions 

<*  © 

g+  and  g  to  obtain 
u ,  v  u ,  v 


and 


E  (r 
z 


a)  g\  , 
y  ,  v* 


0 
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#7b) 


I 

U’.V 


Vr 


a>  V 


0 


(2.48®) 


Equations  2.47  and  2.48  are  zero  over  the  entire  cylindrical  surface  r  * 

a,  a  result  of  g  and  g+  being  defined  only  on  the  helix  tapes, 
°y,v  y ,  v 

while  E_  and  E  are  zero  on  these  same  tapes.  Applying  orthogonality  to 
v 
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Eqs.  2.47  and  2.48  over  the  cylindrical  surface  r  -  a  eliminates  the 

double  summation  over  the  indices  y’  and  v' ;  and  if  the  expressions  in 

Eqs.  2.46a  and  2.46b  are  substituted  for  E_(r  -  a)  and  EJr  -  a),  the 

z  $ 

result  is  four  doubly  infinite  sets  of  linear  homogeneous  simultaneous 

equations  having  the  same  number  of  unknowns.  By  such  manipulations  and 

subsequent  interchanging  the  order  of  the  two  summation  signs  £  and 

a  ,n 

J  ,  one  finally  obtains  from  Eq.  2.47, 
u,v 

■  \ 

l  2j~  i  G;  (ya  aV  kbib  g;  ,  ,  ■ 

L  UtV  0L  £fn;u>v  \l,n  J  l  l  l»n;uvfvf 
y  i  ^  x.  f  n  + 

V  * 

*  l  ZJ*  .Jo*  fir,  a)2  K.I,  of  ,  ,  ■ 

„,v  “•V  l.n  V  l.n;  11 


^  \.v'  l  Gt,n;y,v  ^l6i,na^  KlISL  Gt,n;u*,V 

y ,  v  x.  ,n 


-  I  \J  ■  l  G0  Ubb  a)  K0I0  Gb  ,  ,  -  -  0 

L  y,v  .  i,n;y,v  i,n  t  i  £.,n;y',v* 
y ,  v  ,  n 

(2.49®) 


and  from  Eq.  2.48, 
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A  solution  exists  only  if  the  determinant  of  the  coefficients  of  these 

equations  vanishes.  Thus,  formally  at  least,  a  determinantal  .uation 

can  be  obatained  in  which  8  a  is  calculated  as  a  function  of  ka.  The 

o ,  o 

problem  is  simplified  by  noting  that  the  symmetry  of  the  structure 

♦  - 

allows  J  to  be  determined  from  J  or  vice  versa.  The  four  sets  of 

equations  then  degenerate  into  two  sets  of  independent  equations: 

either  Eqs.  2.*i 9a  and  2.50a  or  Eqs.  2.H9b  and  2.50b. 

However,  the  difficulty  involved  in  the  calculation  of  numerical 

results  is  clear,  and  a  more  convenient  form  for  deriving  approximate 

solutions  is  desired. 

I 
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2.1.6.  Using  the  Variational  Method  to  Derive  an 

Approximation  to  the  Determl nantal  Equation 

As  mentioned  in  Section  1.4.3*  a  variational  technique  Is  useful 
for  deriving  approximate  solutions  to  yield  numerical  results.  Begin¬ 
ning  with  one  of  the  standard  forms  of  the  Lagrangian  for  an  electromag¬ 
netic  field,  a  variational  expression  is  found  for  the  present  problem 
in  terms  of  field  intensities,  which  satisfy  Maxwell’s  equations  as  well 
as  the  symmetry  and  periodicity  conditions,  but  not  the  boundary  condi¬ 
tions  on  the  cylindrical  surface  r  -  a.*  The  result  is  an  expression 
for  the  complex  power,  I,  which  might  be  generated  or  absorbed  by  the 
cylindrical  surface  r  -  a: 
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2m 

dz  J  a  d<p  n 
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(r 
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2H(r 


1 


(2.51  ) 


Because  the  terms  within  the  bracket  are  dotted  with  n,  the  unit  vector 
perpendicular  to  the  surface  r  «  a,  one  needs  only  to  consider  the 
tangential  components  of  E  and  H.  In  light  of  Eqs.  2.16  and  2.17,  Eq . 
2.51  is  rewritten  as 
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*  See  Appendix  C. 
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Writing  the  tangential  components  of  the  H  field  in  terms  of  the  surface 
current  density  (Eqs.  2.2*1  and  2.25)  allows  Eq.  2.53  to  be  expressed  as 
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By  employing  the  usual  Fourier  expansion  for  E  .  E  ,  and  ZJ*  and 
performing  the  integration  (equivalent  to  an  orthogonal ity  integration 
whereby  all  the  cross  terms  of  the  multiplied  summations  are  elimi¬ 
nated),  the  following  form  for  the  variational  expression  results: 


r 

■*  “ 

r 

N 

E»(r  -  s]l,n 

*  ^  ■  <„ 

Zj5  +  Zj0 

l,n  2,,n 

l,n  t,n 

__  mt 

. 

* 

(2.55) 


The  solution  to  the  problem  is  then  given  by 


si  J.  ;  J. 

-!,n  -i,n 


(2.56) 


-  +•  -  + 

but  like  J  and  J  ,  the  small  variations  6J  and  6J  must  themselves 

vanish  off  their  respective  helices.  Once  this  is  taken  care  of,  the 
four  doubly  infinite  set  of  simultaneous  equations  which  result  are 
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Identical  to  Eqs.  2.49  and  2.50. * 

However,  the  goal  here  is  to  find  a  simplified  version  of  the 
deterrainantal  equation.  Using  the  narrow  tape  approximation,  Eq.  2.42, 
the  variational  expression 


I 


is  written  out  in  teri,.o  of  the  E  field  expansions  given  in  Eqs.  2.32  and 
2.33.  After  simplifying  and  making  use  of  the  relations  in  Eqs.  2.44a 
and  2.44b,  this  becomes 


=  [ y  a)2  K  I  sin2e  ♦ 
l,n  V  £,n  I  l  l 


B  s 

r.-H-Jvt  * (ka)" 

t,n 


2. 

cos  8 


-  Uei-na)  KtIt  sin  8 


(2.59) 


f  See  reference  2  for  details,  pp.  47-54. 
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and 


2‘-  ■  M2  v* sln2s  •  (¥) 


2. 

cos  e 


(2.60) 


As  in  Section  2.1.5,  the  surface  current  densities  can  be  made  to 
conform  to  the  constraint  given  in  Eq.  2.36  by  writing  ^  in  terms  of 
a  complete  set  of  orthonormal  functions  which  are  themselves  aefined 
only  on  the  helix  tapes.  However,  for  an  approximate  solution,  may 
be  given  by  a  finite  number  of  terms  of  certain  convenient  functions, 
each  having  an  unknown  coefficient  and  each  defined  only  on  the  helix 
tapes.  Though  the  degree  of  accuracy  increases  as  the  number  of  terms 
in  the  sequence  increases,  it  has  been  found  that  a  one-term  approxima¬ 
tion  provides  good  results,  while  greatly  simplifying  the  variational 
expression.  For  the  one-term  approximation,  there  is  no  variational 
problem,  and  the  determinantal  equation  is  simply 


I  -  0 


Thus,  "j  is  approximated  by 


(2.61  ) 


'  j  g  2 

A  e  0,0  ,  0  S  z  S  p  and  —  -  |  S  <*>  £  —  +  £  (on  helix) 

p  2  P  2 

lj"  -  <  (2.62) 

0  ,  otherwise  (off  helix) 

\ 

The  Fourier  coefficients  are  then  found  by  setting  the  Fourier 

f  n 

expansion  for  lj  (Eq.  2.!J1a)  equal  to  Eq.  2.62  and  apolying  orthogonal¬ 
ity  , 
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(2.63) 


p  -2’2/P‘5/2  eJ‘* 


lj"  m  !  dz  S 

t,n  2mp  q  2irz/p-^/2 


The  integral  over  <(>  reduces  to 


2i rz/p+5/2  Jl 

/  d*  eJU  -  e 

2irz/p-£/2 


2itz 


si 


llii) 


and  the  integral  over  z  then  takes  the  form 


fp  .  J8o.o2  -J6l,n2  Jl  p  e 
Jdze  e  *  e  K  - 

o 


2ttz  Jp(b  -6.  +  )  -  1 

-  JKl  o,o  l,n  p  / 

d  e  n  ' 


O  «  £2ir 

8  -  8,  ♦  - 

oto  I.n  p 


Since 


8.  -  8  +  (l  ♦  2n)  — 

S,,n  o,o  v  '  p 


Eq.  2.65  reduces  to 


"  ■  iik  5l"  <2’n) 


2p  ,  n  -  0 

0  ,  n  *  0 


Considering  Eq.  2.63  in  its  entirety  yields  the  result  that 


,J!,n  -  *  h  *t  5<n) 


where 


(2.6K) 


(2.65) 


2.12') 


(2.66) 


(2.67) 


(2.68) 


140  - 


sin  ( U/2 ) 
(U/2) 


(2.69) 


a 


£ 


and 


6  (n) 


1  ,  n  -  0 
0  ,  n  *  0 


From  Eq.  2.57,  the  determinantal  equation  becomes 


(2.70) 


2U  + 

0,0 


l  al  (v.  •  ■ 0 


(2.71  ) 


The  reason  there  is  no  Z„  term  in  Eq.  2.71  is  because  of  the  delta 

£,n  ^ 


function  in  the  expression  for 


lj- 


i,n 


(Eq.  2.68).  In  the  original 

I 


expression  for  the  variation,  Eq.  2.57,  Z  is  multiplied  by  'J  . 

|  n 

While  n  can  only  be  zero,  the  quantity  (n  +  i.)  is  never  zero  since 
£  Z  1.  This  implies  that 


lj;,nu  ■  A  fe  “1  5("  *  l)  (2'72) 

will  always  be  zero,  as  the  delta  function  will  never  take  the  form 

6(0). 

2.1.7  The  Single  Helix  Determinantal  Equation 

The  corresponding  determinantal  equation  for  the  single  helix  is 
presented  here  for  completeness.  It  is  also  helpful  to  have  it  in  a 
form  which  is  easily  adaptable  to  the  more  general  boundary  value 


problem  described  in  Chapter  Three 


Instead  of  Eq.  2.71,  the  determinantal  equation  for  the  single 


helix  is 


Uo  *  Vo  *  j,  “t  (*l  *  *-l>  ■  0  <2-73> 

Noting  that  Yfl  for  the  contrawound  helix  is  identical  to  Y  for  the 

t ,  O  Jc 

single  helix,  the  terms  UQ,  Y^,  and  Y_^  are  simply  Uo.o-  It,o’  and 
Y  .  ,  respectively.  The  VQ  term  results  from  the  fundamental  component 

“  JC  f  o  ^ 

of  the  TE  fields,  which  is  of  course  not  present  in  the  symmetric  mode 
of  the  contrawound  helix and  has  the  form 

Vq  -  (ka)2  K^(YQa)  r (YQa)  cos20  (2.7^) 

It  should  be  of  no  surprise  that  the  determinantal  equation  for 
the  single  helix  differs  only  slightly  from  that  of  the  twin  helices. 
However,  this  difference  is  enough  to  significantly  alter  the  dispersion 
characteristics,  as  will  be  seen  in  Section  2.2. 


t 


This  property  is  demonstrated  analytically  in  Appendix  D. 


2.2  Results  for  the  Free  Space  Problem 


The  formulation  of  Sections  2.1.6  and  2.1.7  for  the  free  space 

problem  is  implemented  with  a  HP  1000  minicomputer.* 

Fixing  the  values  of  0  and  5,  the  determinantal  equation  is  solved 

numerically  to  obtain  6  a  as  a  function  if  ka,  for  the  contrawound 

O  |0 

helix,  and  6  a  as  a  function  of  ka  for  the  single  helix.  Since  8  a  - 
8Qa,  all  dispersion  plots  are  made  with  respect  to  8Qa.  Furthermore, 
the  ordinate  and  abcissa  are  normalized  in  the  conventional  manner  by 
the  relations 


ka/cot  6  -  p/X 


frees pace 


8a/cot  6  -  p/X 


helix 


where  X  is  the  wavelength. 

Before  proceeding,  it  is  important  to  note  that  the  numerical 
results  compare  well  to  experimental  only  for  narrow  helix  tapes  (£  <  1) 
such  that  the  overlap  region  between  the  two  "touching"  tapes  is  kept  to 
a  minimum.  Nevins  [12]  demonstrates  that  the  discrepancy  between  the 
theoretical  predictions  and  experimental  results  is  due  to  the  currents 
deviating  from  their  respective  helical  paths.  When  current  flow  from, 
one  tape  to  the  other  is  prevented  —  i.e.,  the  two  helix  tapes  are  not 
allowed  to  touch  —  the  predicted  and  experimental  results  compare 
closely  and  are  not  dependent  on  6* 


Refer  to  Appendix  E. 


Figures  2.3  and  2.H  show  how  varying  the  pitch  angle  affects  the 

dispersion,  for  £  -  1  and  £  -  2,  respectively.  The  pitch  angle  is 

reduced  as  cot  6  is  Increased  from  2.5  to  10,  resulting  in  a  decrease  in 

the  group  velocity,  v  (v  -  d(ka)/d(Ba)),  along  with  a  decrease  in  the 

s  g 

phase  velocity,  Vp  (v  -  ka/ga). 

Figures  2.5  and  2.6  show  the  effects  of  varying  £,  while  cot  0  is 
fixed  at  10  and  5,  respectively.  In  both  cases,  changing  the  tape  width 
has  a  negligible  affect  on  the  dispersion  characteristics  for  the  single 
helix,  while  those  for  the  contrawound  structure  are  altered  consider- 
aoly.  The  reason  the  contrawound  circuit  is  so  affected  is  due  to  the 
interaction  between  the  two  helices,  which  becomes  stronger  as  the  tape 
width  is  increased. 

Experimental  results  of  Birdsall  and  Everhart  [U]  are  plotted 
along  with  numerical  results  in  Fig.  2.7.  In  both  cases,  the  dispersion 
increases  with  increased  tape  width.  The  deviation  between  the  theo¬ 
retical  and  experimental  results  is  due  almost  entirely  to  the  fore¬ 
shortened  current  paths  which  result  when  the  two  helix  tapes  are 
allowed  to  touch.  However,  the  effects  of  finite  tape  thickness  must 
also  be  considered. 

Experimental  results  [H]  for  several  tape  thicknesses  are  plotted 
in  Fig.s  2.8  and  2.9.  Correlation  with  theory  improves  as  the  tape 
thickness  is  reduced;  furthermore,  a  comparison  between  the  two  figures 
reveals  that  there  is  better  agreement  between  theory  and  experiment  for 
the  smaller  tape  width,  5  -  tr/H. 


0}O3/e>i 


Fig.  2.3.  Theoretical  dispersion  characteristics  for  both  a  contrawound  helix 
(solid  line)  and  for  a  single  helix  (dashed  line)  in  free  space 
showing  the  effects  of  varying  cot  8:  cot  0  «  2.5,  5,  10.  The 
tape  width  is  described  by  £  =  1  and  is  assumed  to  be  infinitely  thin. 


Fig.  2.4.  Theoretical  dispersion  characteristics  for  a  contrawound  helix 
(solid  line)  and  a  single  heiix  (dashed  line)  of  £  ■  2  in  free 
space,  showing  the  effects  of  varying  cot  0:  cot  9  =  2.5,  5, 
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2.5.  Theoretical  dispersion  characteristics  for  a  contrawound  helix 
(solid  line)  and  a  single  helix  (dashed  line)  of  cot  »  10  In 
free  space,  showing  the  effects  of  varying  C  *  1,  2,  3- 


.25 


Fig.  2.6.  Theoretical  dispersion  characteristics  for  a  contrawound 
(solid  line)  and  a  single  helix  (dashed  line)  of  cot  8  - 
free  space  showing  the  effects  of  varying  £:  £  ■  1 .  2. 
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Fig.  2.9.  Experimentally  derived  dispersion  characteristics  (dot-dashed  line)  for  a  contra- 
wound  helix  in  free  space,  showing  the  effect  of  varying  the  tape  thickness, 
al/ao  =  0.8,  0.9,  are  compared  to  theoretical  results  (solid  line),  in  which 

a(/a  -  1 .  The  other  geometric  parameters  are  cot  8  =  iJ.1*,  £  =  3 it/1*. 


III.  ANALYTICAL  CONSIDERATIONS  OF  METAL  AND  DIELECTRIC 
LOADING  ON  THE  CONTRAWOUND  HELIX 


3.1  The  Boundary  Value  Problem 

To  study  the  effects  of  dielectric  loading  on  the  dispersion 
characteristics  for  the  contrawound  helix,  the  dispersion  equation 
developed  in  Chapter  Two  must  be  altered  so  as  to  allow  for  variations 
in  the  dielectric  properties.  As  before,  the  problem  is  separated  into 
two  regions,  but  each  with  its  own  dielectric  constant.  Figure  3.1 
shows  schematically  that  the  permeability  is  still  that  of  free  space, 
viq,  while  the  permittivity  is  arbitrary,  or  e^.  The  field  quantities 
are  written  as  a  Fourier  decomposition  in  each  region  and  are  then 
matched  across  the  boundary  r  •  a  to  express  them  in  terms  of  the  sur¬ 
face  current  density  on  the  two  helix  tapes.  Finally,  the  variational 
method  is  again  employed  to  obtain  a  determinantal  equation. 

If  Region  2  is  bounded  by  a  conducting  sheath  at  r  -  b,  Fig.  3.2, 
the  problem  becomes  one  in  which  the  contrawound  helix  feels  the  effects 
of  an  external  shield.  This  effect  is  handled  mathematically  simply  by 
reforming  the  modified  cylindrical  Bessel  function(s)  in  Region  2  to 
properly  account  for  this  boundary. 

3.1.1  The  Field  Functions 

Before  writing  the  Fourier  expansions  for  the  field  quantities,  k 
and  y  must  be  redefined  to  correctly  account  for  the  dielectric  proper¬ 
ties  in  each  region.  Thus, 
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erposed  between  two  dielectric  region: 
and  with  permeabilities  u  «  u  -  p  . 
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i  -  1,  2 


(3.D 


ki  ■  “Jeiwo 


(Yn),  '  (6*.r,  -  *fj 


i-  1,  2 


(3.2) 


where  i  denotes  either  Region  1  or  Region  2.  Furthermore,  the  intrinsic 
wave  impedance  is  defined  to  be 


"i-fc  1 '  2 


(3.3) 


To  facilitate  both  clarity  and  understanding,  a  shorthand  notation  is 
adopted  in  which  the  l  and  n  dependence  of  the  terms  in  the  Fourier 
decomposition  for  the  field  quantities  is  assumed: 


'*  ■  (Y..")t 


8*8, 


F(V)  *  r]' 


Fi  ■  h •  Kr  Ki 


(3.M 


(3.5) 


(3.6) 


(3.7) 
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Consequently,  the  Fourier  expansions  may  be  written  in  a  form  which 
assumes  summation  over  the  two  indices  £.  and  n  whereby 


I  ■  I 

l,n--« 


In  this  notation,  the  expressions  for  Ez,  H2,  E^,  and  in  the  two 
regions  are  as  follows:* 


(0  <  r  <  a) 


-JBz  Jl* 


1Ez  -  l  Alt^r)  e  Jte 


1E<t> 


ln.8  k 

—  Al^r)  -  J  ^  Cl'^r) 


Y?r 


e-JBz 


-JBz  U<|> 


H.  '  l 


1  <t> 


niYir 


e-jBz  ejU 


(3.8) 


(3.9) 


(3.10) 


(3.1D 


(a  <  r)  or  (a  <  r  <  b) 


2E«  ‘  1 


2Zz  -  l  n2BKe(V}  eJBZ  eJt* 


ln?B  k 

V-  BK  (Y  r)  -  j  —  DK'(Y  r) 

Y^r  2 


(3.12) 


rJBz  (3.13) 


See  Appenedix  A  for  more  details. 
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-JBZ 


A  - 1  ^  «h<V> e  J  * 


(3. UO 


,H  -  l  -  J  ~  BK*(y2p)  -  -^2“  DKh(Y2r)  •  e~J6z  eJ£*  (3.15) 


2  h  2 
n2T2r 


If  Region  2  is  unbounded,  the  modified  Bessel  functions  in  the  last  four 
equations  are  simply, 


Ke  "  Ki 


[(V")2  '] 


K'  Y.  r 
^  \  ^  * n  /2 


(3.16) 


K'  -  K’  hr.  r 

h  1  \  ♦ny  2 


If,  however.  Region  2  is  bounded  by  a  perfect  conductor  at  r  ■  b,  these 
same  functions  take  the  form  (in  shorthand  notation), 


Kg  -  l(Y2b)  K(Y2r)  -  K(Y2b)  l(Y2r) 


K'e  -  I’(Y2b)  K’(Y2r)  -  K'(Y2b)  I’(Y2r) 


(3.17) 


Kh  -  I’(Y2b)  K(Y2r)  -  K’(Y2b)  I^r) 


-  l(Y2b)  K'(Y2r)  -  K(Y2b)  I’(Y2r) 
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Applying  boundry  conditions  2.1^  and  2.15  at  r  ■  a  and  pr,v'**“di ng 
as  in  Section  2.1.2,  the  Fourier  coefficients  for  the  field  components 
are  determined  in  terms  of  the  surface  current  density  Fourier  coeffici¬ 
ents.  Letting 

zj±  .  ZJ*  (3.18) 

l  ,n 

and 

V  s  (3.19) 

l  ,n 

the  expression  for  C  Is  found  to  be 

C  -  R(ZJ~  +  ZJ+)  +  S(V  +  +  )  (3.20) 


-  58  - 


3 

a  Y  y 
1  2 


wiik;  •  'iva1: 


(3.23) 
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(v)3(v): 


(v)2  •  M 


k2ViiK;  -  vw 


(3.2*0 


and 


-  WlV 


(v)2  -  (v)' 


(3.25) 


Similarly,  the  expression  for  A  is 


A  -  j  £r(zJ*  ♦  V)  ♦  X(V  ♦  V)J  (3.26) 

where 


T  - 


1 

N 


(3.27) 


and 


x--5 


te  . 


M 


181 


(3.28) 


In  the  above  equations,  the  argument  of  I  and  I'  is  Y  a,  while  K^,  K' , 

1  e  e 

K^,  and  take  the  form  of  Eqs.  3.16  and  3.17  with  r  «  a.  Expressions 


for  the  Fourier  coefficients  B  and  D  are  determined  from  the  relation- 


and 


i-  c  -  (V  .  V) 

"l 

3.1.2  Approximating  the  Determinantal  Equation 


(3-30) 


As  outlined  in  Section  2.1.5,  the  solution  to  the  boundary  value 
problem  resulting  in  the  determinantal  equation  can  be  found  by  solving 
four  doubly  infinite  sets  of  linear  homogeneous  simultaneous  equations, 
each  having  the  same  number  of  unknowns.  However,  a  more  manageable 
approach  is  to  use  the  variational  method  of  Section  2.1.6  to  obtain  an 
approximation  to  the  determinantal  equation. 

Since  the  variational  expression  for  the  current  problem  is  the 
same  as  that  for  the  contrawound  helix  in  free  space,  Eq.  2.55  may  be 
used  directly.  From  Eqs.  3.20  and  3.26,  the  expressions  for  E,  and  E 
at  r  -  a  are 


Ez(r  -  a)  -  l  Jn.,1 


T(ZJ~  +  ZJ+)  ♦  X(V  +  *J+) 


e'j6zeJU  (3.3D 


and 
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■JT' 

T(V  *  V)  *  X(V  .  *J*) 

l  V 
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-  j^r 
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r(zj"  ♦  Zj  +  )  ♦  S(V  ♦ 

V  - 

e-JB z  ejU 


(3.32) 
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These  expressions  are  substituted  into  Eq.  2.55  resulting  in  Eq.  2.57, 


with  n,  Y,  ,  and  Z.  defined  as  follows: 
o,n»  £,n  £,n 


U  -  n,IT  sin  e 
o  ,n  1 


Ylfn  -  n^T  Sin^e  - 


tn.B  k 
-J-  IX  ♦  -1  I’S 
T1 


2. 

cos  b 


(3.33) 


and 


tn.  B  k 

n.IX  -  IT  -  — -  I’R 

’  Y^a  1 


sin  e  cos  0  (3.3^) 


Zl,n  ‘  VT  sl"2e  * 


in  6  k 

IX  +  -1  I'S 

Y^a  M 


2 

cos  0 


in.B  k 

n.IX  -  IT  -  —  PR 

'  Y1a  1 


sin  0  cos  0  (3.35) 


Using  the  approximation  for  *J  given  in  Eq.  2.62,  the  determinantal 
equation  is  again  Eq.  2.71,  where  UQ  Q,  Y.  ,  and  Y  are  determined 

'■'t'-'  SC  j  O  “ICjO 

from  Eqs.  3.33  and  3 - 31* . 

Similarly,  for  the  single  helix  interposed  between  two  dielectric 

regions  at  r  >  a  and  with  or  without  a  conducting  sheath  at  r  •  b,  the 

determinantal  equation  is  found  to  be  the  same  as  Eq.  2.73.  As  before, 

UD ,  Y  ,  and  Y  for  the  single  helix  are  simply  U  ,  Y  ,  and  Y  , 

respectively,  but  with  U0  0,  Y  ,  and  Y  being  determined  from  Eqs. 

f  JC  ,  o  -Jc  fo 
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3.33  and  3.3^.  The  corresponding  form  of  VQ  for  the  two  dielectric 
problem  is 


V 

o 


k2n2',i:Ke 


VWh 


e  •  k,"lV 


—  I*  cos  e 


(3.36) 


The  formulation 
explore  the  effects  of 


developed  in  this  section  can 
dielectric  and  metal  loading. 


now  be  used  to 


1 


< 


i 
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3.2  Results  for  Dielectric  and  Metal  Loading 

By  implementing  the  formulation  of  Section  3.1.2,  the  effects  of 
loading  are  investigated  simultaneously  for  both  the  contrawound  and 
single  helix. 

3.2.1  Dielectric  Loading 

The  dispersion  plot  in  Fig.  3-3  for  cot  6-10  shows  how  the  group 
velocity  and  phase  velocity  are  reduced  as  e2/ei  is  increased  from  1  to 
9.  This  behavior  is  also  present  in  the  experimental  work  performed  by 
Birdsall  and  Everhart  [i»].  They  observe  that  as  the  distance  between  a 
surrounding  glass  cylinder  and  the  slow-wave  structure  is  reduced,  the 
degree  to  which  the  velocities  are  loaded  increases.  Their  results  are 
reproduced  in  Fig.  3.^.  Inspection  of  this  figure  further  reveals  that 
the  effects  of  dielectric  loading  are  negligible  when  the  distance 
between-  the  cylinder  and  the  circuit  is  increased  past  a  certain 
point.  The  reason  for  this  is  the  fields  outside  the  circuit  (Region  2) 
decrease  exponentially  and  consequently  do  not  penetrate  radially  a 
significant  distance  into  Region  2. 

Figure  3.5  attempts  to  correlate  the  experimental  results  for 
loading  by  a  glass  cylinder  of  finite  thickness  to  theoretical  predic¬ 
tions  for  dielectric  loading,  in  which  the  dielectric  of  Region  2 
extends  to  infinity.  To  do  so,  an  equivalent  relative  permittivity  for 
Region  2  is  calculated  by  volumetrically  proportioning  the  relative 
permittivity  of  the  glass  cylinder. 

It  is  assumed  that  the  fields  in  Region  2  penetrate  to  a  depth  of 
b/a»  1.21.  Consequently,  the  slow-wave  circuit  "sees"  an  effective  (two 
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Fig.  3-3-  Theoretical  dispersion  characteristics  for  a  contrawound  helix  (solid  line)  and 
a  single  helix  (dashed  line)  interposed  between  two  dielectric  regions  and  show¬ 
ing  the  effect  of  dielectric  loading  by  the  outer  region  --  region  2  --  such  that 
£_/€.  -  1,  1.5,  9  and  k?/k,  «  1,  1.225,  3,  correspondingly,  cot  0  *=  10,  £  «  1 . 


dimensional)  volume  of 


vol  -  — r  (b2  -  a2]  -  t[l .21*  -  ld) 
a 


2  ,2- 


vol  ..  -  1.46 
eff 


Comparisons  can  best  be  made  with  the  experimental  results  (Fig. 

3.4)  for  the  glass  cylinder  with  dimensions  b^/a0  -  1.05  and 

(b  -  b  1/b,  -  0.122.  The  ratio  b„/a^  is  then  1.178  and  the  (two  dimen- 
o  i  i  o  o 

sional)  volume  for  the  cylinder  is  thus 


vol 


gc 


~2  (»o-  »?)-  ’(l-178?  *  1'053) 


vol 


gc 


0.896 


The  relative  permittivity  of  the  low  loss  glass  cylinder  is 
assumed  to  be  3.0,  and  the  equivalent  relative  permittivity  for  Region  2 
is  proportioned  as 


'vol 


3.0 


gc 


vol 


3.0 


eff 


0,896 

1.46 


\ 

V 


1 .83 


The  theoretical  results  for  e_,  =*  1.83  along  with  the  correspond!  ng 

t 

experimental  results  are  snown  in  Fig.  3.5.  The  deviation  between 
theory  and  experiment  is  due  to  the  for eshortened  current  paths  wrich 
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result  when  the  two  helix  tapes  are  allowed  to  touch.  However,  it  is 
clear  that  the  velocity  loading  is  proportionally  similar  in  both  cases. 

Figures  3.6  and  3.7,  like  Fig.  3.3.  allow  e^/e^  t0  vapy  1  to 
9,  while  cot  6  is  fixed  at  5  (Fig.  3.6)  and  2.5  (Fig.  3.7).  In  all 
three  figures,  the  percentage  of  dielectric  loading  relative  to  no 
loading  is  approximately  the  same. 

In  Figs.  3.8  and  3.9,  is  fixed  at  1*5  and  9,  respectively. 
The  overall  shape  of  the  dispersion  curves  changes  as  cot  6  is  increased 
from  2.5  to  10. 

3.2.2  Metal  Loading 

The  effects  of  metal  loading  are  considerably  different  than  the 
effects  of  dielectric  loading.  No  longer  are  there  the  forbidden 
regions  associated  with  the  open  structure.  And  the  general  shape  of 
the  oj-8  diagram  changes  as  the  effects  of  metal  loading  become  stronger. 

The  dispersion  characteristics  for  a  contrawound  helix  (5  -  1, 
cot  6  *  10)  symmetrically  oriented  inside  a  cylindrical  conducting 
sheath  are  presented  in  Fig.  3. 10.  The  relative  constitutive  properties 
are  k-^/k^  -  1,  -  1,  and  e^/e.^  “  1.  By  varying  b/a,  the  relative 
distance  between  the  circuit  and  the  cylinder,  the  aspects  of  metal 
loading  mentioned  above  are  clearly  seen. 

The  dashed  curve,  labeled  1,  is  the  dispersion  for  the  nonshielded 
contrawound  helix.  The  ends  of  this  curve  couple  into  the  so-called 
"velocity  of  light  lines,"  as  this  is  an  open  structure.  Curve  1A  shows 
the  effect  of  a  conducting  shield  placed  radially  at  a  distance  of  b/a  = 
2.  Instead  of  coupling  into  the  velocity  of  light  line,  curve  1A 
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2 25,  3.  corres 
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Theoretical  dispersion  characteristics  for  a  contrawound  helix  of  t;  -  1,  Inter¬ 
posed  between  two  dielectric  regions  with  relative  constitutive  properties  of 
e  /e  =9  and  k^/k,  =  3,  showing  the  effect  of  varying  cot  0:  cot  0  -  2.5,  5, 
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deviates  from  the  nonshielded  case  at  approximately  Boa/cot  6  -  0.8^  and 
couples  into  the  coaxial  TE^  mode.  This  coaxial  mode  is  formed  by  the 
contrawound  helix  as  the  inner  conductor,  and  the  outer  cylinder  as  the 
outer  conductor.  As  the  outer  cylinder  is  brought  closer  to  the  cir¬ 
cuit,  the  cutoff  for  the  coaxial  TE^  mode  moves  up  in  frequency. 
Consequently .coupling  between  the  fundamental  component  of  the  space 
harmonics  and  the  coaxial  TE^  mode  occurs  at  higher  values  of  ka 
(curves  IB  and  1C).  A  point  is  reached  where  the  shield  is  close  enough 
to  the  contrawound  helix  that  coupling  no  longer  occurs  --  the  TE^  mode 
has  moved  out  of  range  (curve  ID). 

As  b/a  continues  to  approach  1,  more  of  the  E-field  between  the 
circuit  and  the  cylinder  is  terminated  on  the  cylinder.  By  this 
mechanism,  the  field  shape  for  the  traveling  wave  remains  unchanged  over 
a  wider  range  of  frequency,  thereby  reducing  the  dispersion  (curves  IE 
and  IF).  For  b/a  «  1,  curve  1G,  all  of  the  E-field  terminates  on  the 
cylinder,  resulting  in  two  dispersionless  helically  traveling  waves  — 
one  right  handed,  the  other  left  handed  —  each  propagating  axially  at 
p/2ira  times  the  velocity  of  light.  Whether  or  not  this  situation  is 
physically  possible  is  uiscussed  in  Chapter  Four. 

As  a  comparison,  the  dispersion  for  the  single  helix  in  free  space 
(nonshielded)  is  plotted  in  Fig.  3.11,  which  repeats  Fig.  3.10  with  this 
addition.  The  dot-dash  single  helix  curve  (labeled  2)  parallels  closely 
curve  1G.  The  deviation  between  the  two  is  shown  on  an  enlarged  scale 
in  Figs.  3.12  and  3.13. 

The  effects  of  metal  loading  on  the  single  helix  with  the  same 
parameters  as  above  are  investigated  in  Fig.  3-11*.  Like  the  contrawound 
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Fig.  3 .12.  Enlarged  portion  of  Fig.  3.11.  Ref 
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Rnlarged  portion  of  Fig.  3.11.  Refer  to  Figs.  3.10  and  3.11  for  detai 
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Theoretical  dispersion  characteristics  for  a  single  helix  of  cot  0-10 
and  £  -  1 ,  symmetrically  oriented  inside  a  conducting  cylinder,  showing 
the  effect  of  varying  b'a.  Dashed  curve  2:  b/a  =  °»  (free  space);  curve 
2A:  b/a  =  2;  curve  2B:  b/a  =  1.5;  curve  2C:  b/a  -  1.1;  curve  2D:  b/a 
=  1.  The  relative  constitutive  properties  are  e  /e  =  1  and  k0/k.  *  1. 


helix,  coupling  with  the  coaxial  TE^  mode  is  reduced  as  b/a  approaches 
1.  And  like  before,  the  ideal  case  of  no  dispersion  is  achieved  when 
b/a  -  1  (curve  2D).  Figure  3.15  shows  how  the  tendency  to  couple  to  the 
velocity  of  light  line  is  reduced  as  the  distance  between  the  shield  and 
the  helix  becomes  smaller. 

The  effects  of  metal  loading  on  a  contrawound  helix  of  cot  0-2.5 
are  presented  in  Figs.  3.16  and  3.17. 

Experimental  analysis  performed  by  Birdsall  and  Everhart  [i»]  on  a 
ring-bar  circuit  inside  a  metal  cylinder  lends  validity  to  the  previous 
theoretical  results.  Experimentally  derived  dispersion  curves  for  two 
different  size  cylinders,  b/a  -  1.33  and  b/a  -  2.16,  are  reproduced  in 
Fig.  3*18.  Also  plotted  are  the  theoretical  predictions  for  the  same 
dimensions.  It  is  clear  that  as  b/a  approaches  1,  the  U-B  curves  become 
less  dispersive.  Observe  also  that  for  values  of  b/a  >  2,  the  effects 
of  metal  loading  are  minimal.  Note  again  that  the  deviation  between 
theory  and  experiment  is  a  consequence  of  foreshortening  the  current 
paths.  For  the  ring-bar  geometry,  this  is  equivalent  to  a  large 
connecting-bar  width. 

3-2.3  Simultaneous  Metal  and  Dielectric  Loading 

How  the  dispersion  is  affected  by  the  simultaneous  loading  result¬ 
ing  from  a  dielectric  region  interposed  between  the  slow-wave  structure 
and  a  metal  shield  can  be  qualitatively  determined  by  simply  superposing 
the  individual  effects  of  dielectric  loading  and  metal  loading.  How- 
I  ever,  consideration  of  accuracy  requires  the  use  of  the  quantitative 

analysis  employed  thus  far. 
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p  1  pr .  3.15.  Fnlarged  portion  of  Fig.  3.1^.  Refer  to  Fig. 
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Theoretical  dispersion  characteristics  for  a  contrawound  helix  of  cot  0-  2.5 
and  £  «  1,  symmetrically  oriented  inside  a  conducting  cylinder,  showing  the 
effect  of  varying  b/a.  Dashed  curve  1:  b/a  =  «*  (free  space);  curve  1A:  b/a 
»  2;  curve  IB:  b/a  -  1.5;  curve  1C:  b/a  «  1.15;  curve  ID:  b/a  -  1.1;  curve  IE: 
b/a  =  1.  The  relative  constitutive  properties  are  e0/e,  -  1  and  k?/k,  -  1. 


3.18.  Comparison  between  theoretically  (solid  line)  and  experimentally  (dot-dashed  line) 
derived  dispersion  characterise cs  Tor  a  contrawound  helix  of  cot  e  -  4.4  and  £  - 
ir/2,  symmetrically  oriented  Inside  a  conducting  cylinder,  showing  the  effect  of  metal 
loading:  b/aQ  =  ®,  1.33,  2.16.  The  relative  constitutive  properties  e  /e  "  1  an d 
k?/k1  =1.  In  the  experiment  aj/a0  -  0.8,  while  aj/a  *  1  is  assumed  in  the  theory. 


In  Figs.  3.19  and  3.20,  the  relative  constitutive  properties  are 


k2/kl  "  1,225»  u2/wl  "  and  c2/cl  "  ^*5*  The  <*»“B  curves  for  various 
values  of  b/a  are  plotted  for  the  contrawound  helix  in  Fig.  3.19,  and 
for  the  single  helix  in  Fig.  3.20.  In  each  case,  the  w-8  curves  become 
ideally  nondispersive  as  b/a  approaches  1.  Concurrently,  the  effects  of 
dielectric  loading  are  reduced  as  the  size  of  Region  2  is  decreased. 
This  behavior  becomes  more  pronounced  as  e^/e^  ^ecomes  larger. 

In  Fig.  3.21,  the  relative  permittivity  is  increased  such  that 
e2/el  *  ®  and  k2/kl  "  3*  The  dispersion  curves  ere  again  plotted  for 
varying  values  of  b/a.  These  plots  are  repeated  in  Fig.  3.22,  with 
the  oj-8  curves  for  a  single  helix  ( £2/el  ’  ^  and  E2/e1  *  ^  added  for 
reference. 
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Fig.  3.20.  Theoretical  dispersici  characteristics  for  a  single  helix  of  cot  8-10  and  5  -  1, 

interposed  between  two  dielectric  regions  symmetrically  oriented  inside  a  conducting 
cylinder,  showing  the  effects  of  varying  b/a.  The  relative  constitutive  properties 
are  e^/e^  -  1.5  and  kp/kj  »  1  .225.  Dashed  curve  2:  b/a  -  »;  curve  2A:  b/a  -  2; 
curve  2B:  b/a  =  1.5;  curve  2C:  b/a  -  1.1;  curve  2D:  b/a  -  1. 00008;  curve  2E:  b/a  -  1 
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IV.  SUMMARY  AND  CONCLUSIONS 


The  formulation  of  Chodorow  and  Chu  [2]  for  the  contrawound  helix 
la  expanded  to  include  the  effects  of  dielectric  and  metal  loading. 
Since  wave  velocity  is  inversely  proportional  to  the  square  root  of  the 
dielectric  constant  whenever  a  dielectric  material  is  added  in  the 
region  surrounding  the  circuit,  the  wave  velocity  is  decreased.  In 
general,  while  dielectric  loading  reduces  both  the  phase  velocity  and 
the  group  velocity,  the  overall  shape  of  the  urB  diagram  remains 
unchanged.  However,  this  is  not  the  case  for  metal  loading. 

As  the  degree  of  metal  loading  is  increased,  the  dispersion  is 
effectively  reduced  for  both  the  single  helix  as  well  as  the  contrawound 
helix.  The  mechanism  whereby  this  is  accomplished  is  provided  for  by 
the  outer  conducting  cylinder.  This  cylinder  allows  the  propagating 
mode  to  retain  In  detail  its  particular  shape  by  providing  an  alterna¬ 
tive  termination  for  the  electric  fields.  Consequently,  the  field 
pattern  for  this  mode  tends  to  change  only  in  scale  as  the  frequency  is 
varied.  This  effect  is  increased  as  the  cylinder  is  brought  closer  to 
the  circuit,  thereby  reducing  the  dispersion.  Unfortunately,  the 
reduced  dispersion  is  offset  by  an  accompanying  decrease  in  the  circuit 
interaction  impedance.  An  inductive  coupling  of  the  helix  currents  to 
the  metal  cylinder  and  the  flow  of  current  in  the  circumferential  direc¬ 
tion  in  the  cylinder  results  in  an  increase  of  the  excess  stored  energy 
in  Er  between  the  circuit  and  the  cylinder,  energy  which  is  useless  for 
Interaction  with  electron  beams. 
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In  the  limit,  the  effects  of  metal  loading  would  be  greatest  if 
the  slow-wave  structure  and  the  outer  cylinder  were  allowed  to  touch. 
The  theory  would  not  break  down  in  this  situation  provided  the  currents 
could  be  maintained  along  their  helical  path.  If  this  could  be  accom¬ 
plished,  the  field  shape  for  the  slow-wave  mode  would  be  "perfectly" 
maintained  independent  of  the  frequency,  and  the  phase  velocity  and 
group  velocity  would  become  p/2ira  times  the  velocity  of  light  —  ideally 
nondispersive.  Unfortunately,  if  contact  were  to  be  made  between  the 
circuit  and  the  cylinder,  the  slow-wave  mode  would  be  shorted  out  such 
that  only  circular  cylindrical  waveguide  modes  could  exist. 
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APPENDIX  A 


THE  FOURIER  DECOMPOSITION  OF  THE  E  AND  H  FIELDS 


For  a  homogeneous  source-free  region,  the  vector  nature  of  an 
arbitrary  electromagnetic  field  can  be  expressed  as  the  sum  of  TE  and  TM 
fields,  and  in  the  case  of  circular  cylindrical  geometry,  these  may  be 
defined  as  being  transverse  with  respect  to  the  axial  coordinate.  For 
regular  boundaries  such  as  that  for  a  circular  waveguide,  the  TE  and  TM 
portions  of  the  fields  uncouple,  giving  separate  solutions.  However, 
for  open  structures  such  as  the  single  helix,  contrawound  helix,  and 
ring-bar  circuits,  the  skew  boundary  conditions  necessitate  both  the  TE 
and  TM  aspects  of  the  fields  to  be  simultaneously  present. 

The  E  and  H,  field  expressions  for  the  free  3pace  problem  shown 
z  z 

in  Fig.  2.1  are 
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and  from  these  expressions  and  Maxwell's  equations,  the  other  field 
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components  take  the  form  in  Region  1 : 
(0  <  r  <  a) 
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and  in  Region  2 
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In  the  above  equations,  the  four  sets  of  Fourier  coefficients 

B„  ,  C.  ,  and  D„  are  the  result  of  the  TE  and  TM  fields  in  the  two 
t,n'  i,n  i,n 

regions,  as  summarized  in  Table  A.1. 

By  expressing  the  Fourier  coefficients  for  the  fields  in  terms  of 
the  Fourier  coefficients  for  the  surface  current  densities  on  the  two 
helix  tapes,  the  field  expansions  for  the  contrawound  helix  take  the 
form: 
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Considered  next  is  the  situation  in  which  Regions  1  and  2  are 
allowed  to  have  individually  distinct  dielectric  properties,  as  illus¬ 
trated  in  Figs.  3.1  or  3.2.  The  intrinsic  wave  impedance  is  given  by 


'■•15 


i  -  1,  2 


(A. 25) 


where  i  denotes  either  Region  1  or  Region  2.  Similarly,  the  radial 
propagation  constant  is  defined  to  be 


(Vn^  *  (8t.n  -  kl)'/2  ‘  ‘  2  <#-26) 

For  convenience,  a  shorthand  notation  is  adopted  whereby  the  l  and  n 
dependence  is  assumed,  In  this  notation, 
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Keeping  in  mind  that  the  Fourier  expansion  of  the  E  and  H  field 

CD 

components  is  over  the  double  summation  £  ,  the  shorthand  form  of 

i , n--® 

these  field  components  becomes  in  Region  1, 
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The  necessary  and  sufficient  condition  for  this  equation  to  be  satisfied 
is 
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where  u  is  an  integer.  Application  of  Eq.  B.17  to  Eq.  B.21  Implies  that 

£  -  n  -  2u  (B.22) 


or 


n  -  2u  +  £ 


(B.23) 


At  this  point,  one  may  pursue  the  analysis  from  two  different  perspec¬ 
tives.  The  first  considers  the  problem  redefined  in  terms  of  a  pseudo¬ 
period  L  -  p/2.  The  double  summation  over  Z  and  n  is  restricted  as  a 
consequence  of  Eq.  B.22  so  that  if  l  is  even,  n  is  even,  and  if  l  is 
odd,  n  is  odd.  This  then  implies  a  coupling  between  $  and  z,  as  would 
be  expected,  and  which  is  necessary  for  the  orthogonality  integral  to  be 
evaluated  over  p/2  rather  than  the  usual  full  period  p. 

The  second  and  more  familiar  scenario  maintains  orthogonality  over 

the  entire  period  p,  but  redefines  the  propagation  characteristics  for 

the  space  harmonics.  Substituting  Eq.  B.23  for  n  in  the  expression  for 

8  given  in  Eq.  B.17  and  noting  u  is  allowed  to  vary  over  the  entire 
n 
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(B.17) 


The  result  is  that  the  propagation  constant  for  the  nth  order  space 
harmonic  is  defined  by  Eq.  B.17. 

When  step-turn  periodicity  (also  called  screw  symmetry)  is  pres¬ 
ent,  as  is  the  case  for  the  contrawound  helix,  Floquet's  theorem  pro¬ 
vides  the  expression, 
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Like  Eq.  B.ll,  the  above  equation  states  that  the  electric  field  evalu¬ 
ated  at  <p  and  z  (written  E($,  z))  and  propagated  a  distance 
-JB  p/2 

p/2  (written  e  )  is  identical  to  the  electric  field  evaluated  at 

the  position  <p  -  0  +  it  and  z  -  z  +  p/2.  The  fact  that  the  fields  at 
these  two  positions  are  inverted  from  one  another  is  automatically 
accounted  for  by  <p  -  $  +  ir. 

By  writing  the  Fourier  expansion  for  the  z  component  of  the  E 
field  at  the  position  $  -  $  ♦  ir  and  z  -  z  +  p/2, 
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and  substituting  this  along  with  Eq.  B.12  into  Eq.  B.18  at  r  -  a,  one 
obtains  the  following  equation: 
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(written  e  0  ),  is  identical  to  the  same  functional  form  for  the  E 
fields  evaluated  instead  at  z  •  z  ♦  p  (written  E(z  +  p)). 

For  a  periodic  slow  wave  structure  with  circular  cylindrical 
geometry  of  radius  "a,"  the  electric  field  can  be  decomposed  into  a 
Fourier  series.  Consequently,  E(z)  is  written  as 
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where  A  is  the  amplitude  factor  and 
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The  Fourier  decomposition  of  E  (r,  z  ♦  p)  is  similarly  given  by 
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Letting  r  ■  a  and  substituting  Eqs.  B.12  and  B.l1!  into  Eq.  B.ll  yields 
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For  Eq.  B.15  to  be  true  implies  that 
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further  implying  that 
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For  a  given  mode  of  propagation  and  at  a  given  frequency, 
the  wave  functions  at  two  points  on  a  transmission  system, 
separated  by  one  period,  differ  by  a  complex  constant. 


The  application  of  Floquet's  theorem  to  periodic  slow  wave  structures 
allows  one  to  determine  the  propagation  characteristic  of  the  space 
harmonics. 

For  a  lossless  periodic  transmission  system,  the  E  fields  are 
giver,  the  form, 
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which  corresponds  to  Eq.  B.7  with  r 
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From  Floquet's  theorem,  it  is  necessary  that  E(z  ♦  p)  ■  E(z),  as  seen  in 
Eq.  B.6.  Thus,  Eq.  B.9  can  be  expressed  as 


"JB  z  -JB  p 

E(z  ♦  p)  -  E(z)  e  °  e 


(B. 10) 
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The  above  equation  simply  states  that  the  functional  form  for  the  E 
fields,  evaluated  at  z  (written  E ( z ) )  and  propagated  a  distance  p 
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A  useful  corollary  to  Floquet's  theorem  is  obtained  in  the  follow¬ 


ing  manner.  Let  k  take  the  form, 

,  npr 
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If  *(x)  is  defined  as 
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then 
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Applying  Eq.  B.2  to  the  above  equation  gives 
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so  that  Floquet's  theorem  states  that  one  can  always  find  a  solution  to 
Mathieu's  equation  of  the  form, 


y(x)  *  e^x  $(x) 


(B.  7) 


where  *(x)  is  periodic  with  period  p. 

It  should  be  apparent  from  Eq.  B.7  that  the  functions  y(x)  can 
represent  wave  functions.  In  light  of  this,  Floquet's  theorem  becomes  a 
statement  about  periodic  translational  symmetry.  For  a  periodic  trans¬ 
mission  system,  Floquet’s  theorem  may  be  stated  as  follows: 
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APPENDIX  B 


FLOQUET’S  THEOREM  AND  STEP-TURN  PERIODICITY 

Floquet’s  theorem  results  from  a  consideration  of  the  second  order 
linear  differential  equation, 


where  a  and  b  are  real  constants.  Known  as  Mathieu's  equation,  Eq.  B.l 
occurs  in  problems  of  wave  motion  with  elliptical  boundaries,  the 
simplest  example  being  the  vibrations  of  an  elliptical  drum  head. 
Although  the  equation  can  be  solved  by  the  usual  power  series  method 
(method  of  Frobenius),  such  a  solution  is  not  valid  when  x  •  np,  "n" 
being  an  integer.  Other  methods  are  then  employed  to  arrive  at  the  more 
general  solution  and,  of  particular  interest,  solutions  which  are 
periodic  in  x,  i.e.,  y(x)  -  y(x  +  np),  where  n  is  an  integer  and  p  is 
the  period.  However,  such  periodicity  is  obtained  only  when  the  con¬ 
stant  "a"  is  allowed  certain  values.  If  "a"  differs  from  these  allowed 
values,  then  the  solution  is  no  longer  periodic.  Instead,  it  takes  the 
form 


y(x  ♦  np)  -  ky(x)  (B.2) 

k  being  a  (complex)  constant.  Equation  B.2  is  the  statement  of  Flo- 
quet's  theorem  and,  though  not  performed  here,  it  is  a  simple  matter  to 
prove  this  theorem  using  linear  algebra  [6,  7]. 
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Finally,  the  Fourier  coefficients  B  and  D  can  be  ♦‘/pressed  in 
terms  of  and  zj±  through  the  relationships 
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At  r  *  a,  E  and  E  in  terms  of  the  surface  current  density  are 
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If  n1  *  n2«  Eqs.  a. 60  and  A. 61  reduce  to  Eqs.  A. 13  (A. 14)  and  A. 15 


(A. 16)  for  r  -  a,  respectively. 
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A  conducting  sheath  placed  at  r  -  b  in  Region  2  (Fig.  3.2)  effects  only 
the  last  four  equations.  These  become  (in  shorthand  notation) 
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The  expression  for  A  is 
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In  the  above  equations,  the  argument  of  the  modified  Bessel  functions  I 
and  I*  is  a,  i.e., 
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where  K  ,  K' ,  Kh,  and  K'  have  the  functional  form  given  in  either  Eqs. 

e  e  n  n 

3.16  or  Eq.  3.17.  As  before,  the  Fourier  coefficients  A,  B,  C,  and  D 
are  a  consequence  of  the  TE  and  TM  fields,  as  outlined  in  Table  A.l. 

Next,  the  Fourier  coefficients  A  and  C  can  be  written  in  terms  of 
the  z  and  $  components  of  the  surface  current  density  on  the  two 
helices.  Letting 
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and 


l  ,n 


(A.HH) 


the  expression  for  C  is 


C  -  R(V  +  V)  ♦  S(V  ♦  V) 


(A.  145) 


where 


R.-ll 

R  n2  IN 


(v)'  -  M 


(A. 146) 


~  Q_  M. 

1  “  *2  IN  Y?a 


1  u 


(v)2  •  (v); 


k1 1  'Yj3  2aj  ni 

3 


£8n2I- 


2M 

N 


M2  -  (v): 


♦  Q 


(A.H7) 
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since  E  ■  0  at  r  »  •»,  and  from  No.  2, 


/  E  •  s  •  E  ds  »  J  E  •  n  »  V  x  E  ds  ■  0  (C.9) 

helix  tapes  helix  tapes 

A 

since  nxVxE-Oon  tne  helix  tapes.  Similarly,  applying  No.  3  to 
the  infinite  planes  at  z  -  0  and  z  -  p  gives 


*  •  2tr 

J  E  •  s  •  E  ds  -  J  r  dr  J  d$ 

planes:  o  o 

z-0,  z-p 

«  ~-0 
♦  e  E  (z  -  0)  •  s(z  -  p)  •  e  E(z  -  0) 


E  (z  ■  0)  •  s(z  -  0)  •  E(z  *  0) 
“JB.P 


“  2lT  # 

«  J  rdrj  d#  E  (z  -  0) 

o  o 


s(z  -  0)  +  s(z  -  p) 


E(z  -  0)  -  0 

(C.10) 


since  s(z  -  p)  -  -  s(z  -  0)  on  account  of  a  reversal  in  the  direction 

A 

of  n.  Considered  together,  Eqs.  C.8  through  C.10  imply  that 


j  E*  •  s  •  E  ds  -  0  (C.ll) 

S 


and  Eq.  C.l  reduces  to 


L 


(V 


,  2  * 
f)  -  k  E 


dv 


(C . 12) 


By  application  of  the  vector  Green’s  theorem  to  Eq.  C.12,  L  is  trans¬ 
formed  to 
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to  be  satisfied  requires  that 


7  x  7  *  E  -  k2E  in  V  (C.6) 

A 

n  »  7  «  E  *  s  •  E  onS  (7.7) 

Thus,  E  must  satisfy  the  vector  wave  equation  in  Eq.  C.6  and  the  bound¬ 
ary  condition  in  Eq.  C.7,  as  expected. 

Next,  for  the  boundary  value  problem  consisting  of  the  contrawound 
helix  in  free  space,  V  is  taken  to  be  an  infinitely  large  circular 
cylinder  having  infinite  radius  and  finite  length  of  one  helix  pitch 
p.  The  surface  S  enclosing  V  consists  of  the  two  infinite  planes  z  «  0 

and  z  -  p;  one  cylindrical  surface  at  r  ■  »,  0  i  z  S  p;  and  the  metallic 

surface  of  the  helix  tapes  at  r  -  a. 

The  physical  nature  of  the  problem  is  such  that  E  must: 

1.  Vanish  at  r  - 

2.  Have  no  tangential  component  on  the  metallic  surface  of  the 
helix  tapes. 

3.  Satisfy  Floquet's  for  a  periodic  structure,  i.e., 

-JB  p 

E(z  -  p)  -  e  °  E(z  «  0) 

Using  these  conditions,  the  surface  integral  in  Eq.  C.l  can  be  shown  to 
be  zero.  From  No.  1, 


J  E*  •  s  •  E  ds  -  0  (C.8) 
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6L  -  /  (7«  6E )*  •  (7  x  E)  dv  ♦  J  (7  x  E)*  •  (7  «  fiE)  dv 

V  V 


2  r  *  2  f  * 

-  k  J  6E  •  E  dv  -  k  J  E  •  5E  dv 

V  ’  '  V 


J  6E  •  s  •  E  ds  +  J  E  •  s  •  5E  ds 
S  S  ‘ 


(C.2) 


or  since 


J  E  •  s  •  5E  ds  -  J  6E  •  s  •  E  ds 
S  '  ~  S 


(C.3) 


Eq.  C.2  is  written  as 


6L  ■  J  (7  x  6E)  •  (7  x  E)  -  k^6E  •  E  dv  ♦  j  6E  •  s  •  E  ds  +  c.c. 

V  L  ‘  “J  3 


(C.4) 


Note  that  the  unexpressed  terms  in  Eq.  C.iJ  are  extraneous  for  the  fol¬ 
lowing  discussion. 

By  means  of  the  vector  Green's  theorem,  6L  can  be  transformed  to 
the  form, 


6L  -  /  6E  •  (7  x  7  x  E  -  k2E)  dv 
V 


-  j  6E  •(nx7xE-s,E)ds+c.c.  (C.5) 

S 


Because  the  first  variation  of  L  is  zero  and  6E  is  an  arbitrary  varia¬ 
tion,  each  integral  in  Eq .  C.5  is  individually  zero.  For  this  condition 


-  in  - 


APPENDIX  C 


DERIVATION  OF  THE  LA GRAN GIAN  USED  IN  CONJUNCTION  WITH 
THE  VARIATIONAL  METHOD  OF  SOLUTION' 


To  employ  the  method  of  calculus  of  variations  (the  variational 
method),  it  is  necessary  to  determine  the  correct  form  of  the  Lagrangian 
for  the  problem.  One  begins  with  the  general  form  of  the  Lagrangian  for 
Maxwell's  equations: 


(y  »  e)  •  (y  *  e) 


2  * 
k  E 


dv  ♦ 


J  s’ 

S 


E  ds 


(C.l) 


In  Eq.  C.l,  S  is  the  surface  enclosing  the  volume  V  under  consid¬ 
eration.  The  most  general  boundary  condition  on  S  is  represented  by  s 

a 

such  that  nxV*E-s*E.  Conditions  2.14  and  2.15  (noting 

A  A 

n  x  V  k  e  »  -JcDun  x  H)  are  then  a  specialized  case  of  s  •  E.  It  should 

be  pointed  out  that  s  is  Hermitian  and  that  the  relationshi  s  •  E  « 

* 

E  •  s  is  valid  for  any  arbitarary  vector  E. 

The  Lagrangian  must  be  constrained  to  match  the  boundary  value 
problem.  Consequently,  the  first  variation  of  L, 

6L  -  0 

is  performed  giving 


t 


See  also  reference  2,  pp.  47-51 


(B.29) 


B  ♦  (l  + 
o  p 


3n) 


The  single  helix  is  invariant  under  the  differential  step-turn 
symmetry, 

E(r,  ♦  •  f- 

The  coupling  expression  for  £  and  n  (equivalent  to  Eq.  B.22  in  the  case 
of  180°  step-turn  symmetry)  is  found  to  be 


Az,  z  ♦  Az 


)■ 


-J60dz 


E(r,  <p,  z) 


(B.30) 


£ 


2p 

n  »  ~  u 
Az 


(B.31) 


Though  u,  like  £  and  n,  is  allowed  to  vary  over  the  entire  range  of 
integers,  Eq.  B.31  is  true  for  all  values  of  Az  only  if  u  «  0.  There¬ 
fore,  £  -  n  such  that  the  Fourier  decomposition  of  Ez  for  the  single 
helix  needs  to  be  summed  over  only  one  index. 


I 
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range  of  integers,  one  can  define  0 


to  be 


8,  -  6  *  fp-  a  ♦  2u) 

i,u  o  p 


(B.24) 


or,  since  u  and  n  are  dummy  indices, 


8.  -  8  ♦  —  (i  *  2n) 

l,n  o  p 


(B. 25) 


The  Fourier  expansion  for  E_  is  then  rewritten  as 


_  ,  r  .  MYfc.nr^  “JYnZ  .11* 

E*<r>  *’  ‘  tl  **•"  e 


(B. 26) 


in  which  6.  is  given  by  Eq.  B.25  and  Y  is  defined  by 
5c  f  n  sc  i  n 


2  2  2 
Yn  ‘  Yn  '  k 


( B. 27 ) 


Thus,  coupling  between  <f>  and  z  occurs  through  the  propagation  constant 
for  the  t,nth  order  space  harmonic. 

It  becomes  apparent  that  if  step-turn  symmetry  is  present  in  a 
given  problem,  then  coupling  exists  between  *  and  z  through  the  propaga¬ 
tion  constant  0^  For  example,  the  120°  step-turn  symmetry 


/  n\  ‘JB  P/3 

E  (  r,  <|>  ♦  Y  •  z  *  §)  "  e  E(r*  z) 


(B. 26) 


results  in  0.  being  given  by 
t,n 


-  in  - 


L  -  /  E*  •  (V  x  V  *  E2  -  k2E2)  dv 


♦  J  E*  -  (V  *  V  x  E2  -  k2E2)  dv 


f  * 

-j  E  •  n  *  7  «  E  ds 
S 


t  #  ^  A 

-  j  E  •  (n^  xV*E^+n2*Vx  E^)  ds  (C.13) 

r-a, 

excluding  tapes 

where  the  subscripts  1  and  2  denote  quantities  for  r  S  a  and  r  2  a, 
respectively.  The  first  three  integrals  in  the  above  equation  are  zero 
because  of  Eqs.  C.6  and  C.ll,  and  the  fourth  can  be  extended  over  the 
entire  surface  r  -  a  because  of  Eq.  C.9.  With  these  simplifications, 
the  Lagrangian  is  restricted  to  an  integral  over  the  surface  r-a.  Let 
I  denote  the  Lagrangian,  now  properly  constrained  to  fit  the  boundary 
value  problem.  Thus,  from  Eq.  C.13, 


I 


-  J 

r-a 


* 

E 


(n^  >«  V  x  E^  n^  x  V  *  )  ds 


(C.1H) 


»  Juju  J  E  •  n^  x  (h^  -  ds  (C.15) 

r-a  ~  ^ 

-  '  + 

in  which  the  relation  -JujuH  -  7  *  E  is  employed  and  noting  -  -  n^. 


*  To  facilitate  understanding,  the  problem  has  been  limited  to  the  case 

°f  W1  “  v2  "  v 
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By  expressing  in  terms  of  surface  current  densities,  one  obtains 

the  variational  integral 


I 


r 


- 

1 

a) 

V  *  V 

# 

♦  E 2(r  -  a) 

V  ♦  V 

-  J 

L  J 

* 

(C.16) 


The  exact  solution  to  the  boundary  value  problem  will  have  I  -  0.  This 
is  because  the  integrand  of  I  vanishes  everywhere;  E  (r  -  a)  and 
Ez(r  -  a)  vanish  on  the  helices,  while  and  ZJ±  vanish  off  the 
helices.  However,  the  converse  statement  is  not  true,  and  to  obtain  the 
exact  solution,  one  must  solve  the  variational  equation 

61-0  (C .17) 


Nevertheless,  one  can  get  an  approximate  solution  by  simply  solving  the 
equation  I  -  0. 

The  physical  interpretation  of  I  is  very  simple;  I  is  the  complex 
power  which  might  be  generated  or  absorbed  by  the  cylindrical  surface 
r  -  a.  It  is,  therefore,  certainly  reasonable  that  I  should  be  zero. 

The  case  of  a  contrawound  helix  separating  two  dielectric  regions 
at  r  -  a  and  bounded  by  a  perfectly  conducting  cylinder  at  r  -  b  is 
analyzed  in  much  the  same  fashion.  The  resulting  variational  integral 
is  again  C.16,  which  is  not  surprising  based  on  the  physical  interpreta¬ 
tion  of  I. 
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APPENDIX  D 


RELATIONSHIPS  BETWEEN  THE  COEFFICIENTS  OF  THE  FOURIER 
DECOMPOSITION  OF  THE  SURFACE  CURRENT  DENSITIES 
FOR  THE  SYMMETRICAL  MODE 


The  symmetrical  relationships  between  the  parallel  and  perpendicu¬ 


lar  components  of  the  surface  current  densities  are 


z ) 


(2.40a) 


V(a, 


z)  - 


V(a.  - 


4>»  z) 


(2.40b) 


the  superscripts  and  "  +  "  representing  the  right-handed  and  left- 
handed  helices,  respecti vely .  With  the  Fourier  decomposi tion  of  each 
surface  current  density  component  given  as 


<t>,  z) 


eJH 


(2.41a) 


J  (af  z)  - 


!.,n 


(2.41b) 


one  can  determine  a  relationship  between  the  two  sets  of  coefficients 

I  ♦  I  -  x  * 

■J.  and  'J„  ,  as  well  as  the  two  sets  of  coefficients  J.  and 

i,rt  t,n  £,n 


Ja 

°  t  n 


th 


Since  the  propagation  constant  for  the  t,n  order  space  harmonic 
is  defined  as 
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(2.12) 


8.  -  8  ♦  —  [i  ♦  2n] 

£.  ,n  op1  1 


it  is  observed  here  that 


8  .  „  ■  8  +  —  [-1  +  2(n  ♦  t) ]  -  8. 

o  p  L  *.,n 


( 2. M3a ) 


and  that 


'  B0  [‘  *  ?<"  -  *>)  '  B., 


(2.143b) 


Next,  Eq.  2.H0a  can  be  written  in  terras  of  the  Fourier  expansion  given 
in  Eq.  2.41a: 


I  lj*  e"J6l’"Z  -  I  lj;  e'JSl'"Z 

tin  1>"  t,n  l’n 


(D.l) 


Letting  l  *  -  l  and  n  ■  n  +  l  for  the  right-hand  side  (RHS)  of  Eq.  D.l 


gives 


I  h;  e*J8l'"Z  eJ‘*  -  I  L-.  ,  e'SKl’^lZ  e}l*  W.2) 

t,n  l-n  tin 


Note  that  the  RHS  still  sums  over  the  same  range  as  the  left-hand  side 
(LHS)  as  a  consequence  of  l  and  n  spanning  the  entire  set  of  integers. 
Substituting  8.  for  8.  (Eq.  2.43a)  in  the  above  equation  finally 

1 1  fi  i  % 

gives  the  result, 


I  h:  e'J8*’"Z  eJt«  -  I  lj-  .  e'J9l-"Z 
.  £,n  L  -l,n+l 

l  ,n  ’ 


(D .  3) 
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From  orthogonality,  the  only  way  the  LHS  can  equal  the  RHS  is  if  the 
coefficients  are  equal  term  by  term.  Thus, 


I 


lr 


-£,n+£ 


(2.Wa) 


Referring  again  to  Eq.  D.l,  one  could  just  as  easily  make  the  substitu¬ 
tion  £  *  -  £  and  n  *  n  +  £  for  the  LHS  to  produce  the  relation, 

i 


-£ , n+£ 


Ir 


£  ,n 


(2. 44b) 


Similar  relationships  hold  true  for  the  perpendicular  components. 
Namely,  that 


and 


-£ ,n+£ 


(D.H) 


X  + 


-£,n+£ 


(D.5) 


Furthermore,  using  the  relationships  in  Eqs.  2.38  and  2.39,  it  is  a 
simple  matter  to  show  that 


_  _  *j- 

£,n  -£,n+£ 


(D .  6 ) 


The  expression  for  Hz  as  given  in  Eq.  2.3^  is  rewritten  here  as 
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r  i  a 


H 


*  l,n—  *'n 


:H.  ■  -  (y,  a)  (*J?  ♦  ^ 

l,n  K  l,n  \  i,n  £,n J 


✓  \ 

I I ( Y.  a)  Kfl(YB  r) 

V  i,n  l  i,n 

< 

/ 

►  eJ£* 

V  *  / 

N 

r  S  a 


(D.7) 


(D.8) 


By  substituting  Eq.  D.6  into  Eq.  D.8,  one  finds  that 


SH0  -  -  [ya  a)  (*J~  -  *j"  .)  •  0  for  l 

l  ,n  £,n  ^  4,n  -i,n+£J 


0  (D.9) 


Equation  D.9  shows  that  the  fundamental  component  of  any  wave  field 

operating  in  the  symmetrical  mode  is  a  pure  TM  field;  the  TE  parts  of 

the  fundamental  component  arising  from  the  two  helices  cancel  each 

other.  In  other  words,  the  symmetrical  solution  has  no  H  component 

~z 

with  4-0.  This  is  rigorously  true  because  no  approximation  is 
involved. 


APPENDIX  E 


NUMERICAL  CONSIDERATIONS 

The  numerical  solutions  to  the  determinantal  equations  are 
obtained  with  a  Hewlett  Packard  1000  minicomputer. 

Because  of  numerical  limitations  and  because  the  series  in  Eqs. 
2.71  and  2.73  converge  slowly,  these  equations  are  reformed  to  provide 
rapid  convergence  [8,  9].  Equation  2.71  is  transformed  to 


where 


L  •  (y  a)2  sir.20  -  (ka)2  cos20 
o.o  \o,o  J 


(E.2) 


and 


A  U)  *  I  C0P4^  *  1  •  2002  +  ^  £2  logCO 
3  i«1  l3 


.  1  _  _ JL  u  . 

U  *  288  ** 


A  (o)  ■  l  ~  -  1  .2002 

3  1-1  l3 


(E.3a) 


(E. 3b) 


In  Eq.  E.l,  the  series 
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(E.4) 


m 

l 

1-1 


2  sin  e 

°l  Lo,o  £ 


is  chosen  because  It  converges  at  about  the  same  rate  as  the  series 


(E.5) 


The  term  LQ  Q  is  reasoned  by  examining  the  quantity  (Y^  q  ♦  Y_^  for 
1-0.  From  the  physical  considerations,  the  expression 


2Y 


0,0 


> 


(E.6) 


is  guaranteed  for  all  l. 


Thus,  L.  _  is  a  consequence  of 

O  (V 


L 

0,0 


2Yn  -  2 

0,0 


(y  a] 2  K  fy  .  a^  I  (y  a  sin20] 

\  o,o  /  °V°'0  /  0  \  0,0  J 


(ka)2  K'fy  a)  I'fy  a)  cos2e 

o\o,oJ  0^0,0  j 


(E.7) 


Equation  E.2  results  from  making  the  approximations  that  KqIo  «  1/2  and 

K'  I '  -  -  1/2,  for  Y  a  -  0(1)  . 
o  o  0,0 

2 

The  series  E.M  summed  by  manipulating  a^/i.  into  a  suitaWe  form 
using  the  definition  for  given  in  Eq.  2.69.  This  form  is  as  follows: 


sin2(U/2) 

1 V 


(U/2); 


1 _ cos(lE) 
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or 


•  a2  f  -i 

J,  r  •  p  [a3(0)  ■  V5)J  <e-8) 

where  is  defined  by  Eq.  E.3. 

Using  the  same  reasoning,  Eq.  2.73  is  transformed  to 

U°  *  V°  *  Lo  Sln  e  2~2  [v°>  '  V5)]  *  l  4  (*i  *  -  E0  -  0 

(E.9) 


where 


Lo  *  (Yo3)2  Sin20  "  ^ka)2  cos  6  (E. 10) 

Equations  E.1  (resp.  Eq.  E.9)  is  solved  numerically  to  obtain  ka 
as  a  function  of  Bo^Qa  (resp.  BQa )  for  the  contrawound  helix  (resp. 
single  helix)  in  free  space.  For  the  two  dielectric  problem  of  Chapter 
Three,  LQ>0  (resp.  L  )  is  scaled  by  the  factor 


in  which  the  argument  of  the  modified  Bessel  function  is  y  a  (resp. 

l,o  1  y ’ 

V)- 

To  further  facilitate  numerical  calculation,  a  numerical  device  is 
employed  which  extrapolates  an  infinite  series  to  its  true  sum,  using  a 
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finite  number  of  approximations.  For  example,  if  f1  represents  a  par¬ 
tial  sum  of  Eq.  E.  1 , 


f 


i 


2U 


0,0 


♦  L  sin  6 
0,0 


+ 


i 


I 

£-1 


(E.11) 


then  the  true  value  of  the  infinite  sum  is  approximated  by  the  following 
extrapolation: 


f 


f i+2  ”  /f 


(fl+2  ~  fi+l) 


-  2f 
i+2  i+1 


(E. 12) 


2 

Known  as  Aitken*  s  6  process  [11],  this  extrapolation  may  be  used  to 
accelerate  the  convergence  of  linear  iterations  provided 


(f-  -  fI*l)  ‘  Cl(f-  -  fl)  lCil  <  1  <E-'3> 

The  physical  considerations  again  ensure  that  condition  E . 1 3  is  satis¬ 
fied  for  the  determinantal  equations  E.1  and  E.9. 
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MISSION 

Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (C3I)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C5/  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability  /maintainability  and  compatibility. 


